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Department  of  Defense 
World  Geodetic  System  1984 

Its  Definition  and  Relationships 
with  Local  Geodetic  Systems 

FOREWORD 

1.  This  technical  report  presents  the  Department  of  Defense  (DoD) 
World  Geodetic  System  1984  (WGS  84)  .  The  development  of  WGS  84  was 
initiated  for  the  purpose  of  providing  the  more  accurate  and  updated 
geodetic  and  gravitational  data  required  by  DoD  weapon  and  navigation 
systems.  The  present  WGS  represents  the  Defense  Mapping  Agency  (DMA) 
modeling  of  the  earth  from  geometric/  geodetic,  and  gravitational 
standpoints  using  data,  techniques,  and  technology  available  through 
early  1984,  However,  the  datum  transformation  relationships  with 
geodetic  datums/systems  have  been  updated  and  revised  based  on 
information  available  through  early  1991. 

2.  DMA  TR  8350.2  contains  no  copyrighted  material,  nor  is  a 
copyright  pending.  Distribution  is  unlimited.  Copies  may  be  requested 
from  DMA  as  indicated  in  the  PREFACE. 
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PREFACE 


This  technical  report  presents  the  Department  of  Defense  (DoD) 
World  Geodetic  System  1984  (WGS  84).  The  major  additions  and 
modifications  to  the  second  edition  include  new  transformation 
constants  for  geodetic  datums  and  reference  systems,  deletion  of 
multiple  regression  equations  for  small  and  isolated  areas,  and  changes 
in  symbols  for  the  ellipsoidal  and  orthometric  heights.  In  an  effort  to 
make  this  report  a  complete  entity  on  its  own,  the  most  important  and 
frequently  used  information  from  its  supplements  has  been  merged  and 
included  in  the  second  edition.  Thus,  this  modification  has  made  it 
possible  to  eliminate  the  Supplement  Part  II  (DMA  TR  8350. 2-B)  and, 
henceforth,  there  will  be  only  two  supplements  to  this  report. 

Supplement  Part  I  (DMA  TR  8350. 2-A)  discusses  WGS  84  and  the 
methods,  techniques,  and  data  used  in  developing  the  parameters  and 
products  defining  it.  Considerable  space  is  devoted  in  Part  I  to  the 
discussion  of  the  WGS  84  Reference  Frame,  Ellipsoid,  Ellipsoidal 
Gravity  Formula,  Earth  Gravitational  Model,  Geoid,  and  methods  and 
procedures  for  obtaining  WGS  84  coordinates.  There  is  no  change  to  this 
supplement . 

Supplement  Part  III  (DMA  TR  8350. 2-C)  comprises  the  classified 
information  for  WGS  84.  However,  the  associated  Earth  Gravitational 
Model  coefficients  above  degree  (n)  and  order  (m)  18  and  the 
corresponding  geoid,  previously  classified,  have  now  been  declassified. 
This  supplement,  which  is  still  classified,  will  be  renumbered  as  Part 
II  when  reprinted  in  the  future. 

Also  distributed  with  the  technical  report  is  a  software  program 
for  datum  transformation  and  coordinate  conversions  called  MADTRAN- 
edition  2..  The  MADTRAN  program  (for  Mapping  Datum  Transformation)  is 
provided  on  5.25  inch  floppy  disc  (double  density)  for  IBM  compatible 
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personal  computers.  The  program  allows  input  from  geodetic,  Universal 
Transverse  Mercator  (UTM) ,  or  the  Military  Grid  Reference  System  (MGRS) 
coordinates.  Over  100  datums  are  available  for  transformation  to  or 
from  WGS  84.  Output  is  automatically  presented  as  geodetic,  UTM,  and 
MGRS  coordinates. 

Users  recjuiring  any  specific  information,  or  any  clarification,  or 
data,  should  contact: 


Director 

Defense  Mapping  Agency 
ATTN:  PR,  ST  A-13 
8613  Lee  Highway 
Fairfax,  VA  22031-2137  (USA) 

Similarly,  requesters  requiring  the  positioning  of  sites  of 
interest  directly  in  WGS  84  via  satellite  point  positioning  should 
contact  the  above  address.  Other  WGS  84  related  requests  and/or 
questions  may  also  be  referred  there. 

Since  WGS  84  is  comprised  of  a  consistent  set  of  parameters,  other 
DoD  organizations  should  not  make  a  substitution  for  any  of  the  WGS  84 
related  parameters/equations  in  an  attempt  to  improve  the  accuracy. 
Such  a  substitution  may  lead  to  less  accurate  WGS  84  products  and  may 
have  adverse  effects. 
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Copies*  of  this  technical  report  may  he  requested  from: 

Director 

Defense  Mapping  Agency 
Combat  Support  Center 
ATTN:  PMSR,  ST  D-17 
6001  MacArthur  Boulevard 
Eethesda,  MD  20816-5001  (USA) 

Phone:  (301)  227-2534 
1-800-826-0342 

*  Note;  Non-DoD  users  can  obtain  copies  of  this  report  at  cost.  Call 
the  above  number  for  further  information. 
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1. 


INTRODUCTION 


The  Defense  Mapping  Agency  (DMA)  produces  numerous  mapping, 
charting,  geodetic,  gravimetric,  and  digital  products  in  support  of  the 
Department  of  Defense  (DoD) .  It  is  advantageous  to  refer  these  products 
to  a  single  geocentric  coordinate  system  for  many  reasons  other  than 
ease  of  working  with  a  large  number  and  variety  of  systems .  Such  a 
system  is  needed  due  to  accuracy  and  user  interface  considerations,  the 
need  for  a  product  to  support  the  widest  possible  range  of  applications 
(local,  worldwide) ,  the  need  to  relate  information  from  one  product  to 
data  obtained  from  another  source  (e.g.,  map/chart  positions  to 
coordinates  obtained  from  inertial  navigation  systems  in  real  time) , 
and  the  need  to  ensure  a  smooth  transition  in  product  use  from  one  part 
of  the  world  to  another. 

In  accomplishing  the  preceding,  such  a  geocentric  system,  termed  a 
world  geodetic  system,  provides  the  basic  reference  frame  and  geometric 
figure  for  the  earth,  models  the  earth  gravimetrically,  and  provides 
the  means  for  relating  positions  on  various  geodetic  datums  and 
systems  to  an  Earth-Centered,  Earth-Fixed  (ECEF)  coordinate  system.  In 
brief,  a  world  geodetic  system  serves  as  the  framework  for  DMA  products 
and  worldwide  DoD  operations. 

Previously,  three  such  systems.  World  Geodetic  System  1960  (WGS 
60),  WGS  66,  and  WGS  72,  each  successively  more  accurate,  have 
supported  DoD  activities.  Although  WGS  72  has  aged  gracefully  and  is 
still  adequate  for  some  DoD  applications,  it  has  several  shortcomings 
which  negate  its  continued  use.  For  example,  the  WGS  72  Earth 
Gravitational  Model  (EGM)  and  Geoid  are  obsolete  and  local  geodetic 
datum-to-WGS  datum  shifts  of  improved  accuracy  and  greater  geographic 
coverage  are  needed  than  are  available  from  WGS  72.  In  addition, 
relatively  minor  orientation  and  scale  errors  also  affect  WGS  72.  Other 
factors  contributing  to  the  desirability  of  replacing  WGS  72  with  an 
improved  system  are: 


1-1 


-  Such  an  updai.e  and  replacement  occurs  at  a  time  when  other 
geodetic  system  changes  are  either  underway  or  contemplated; 
e.g.,  the  up-dating,  readjustment,  and  replacement  of  North 
American  Datum  1927  (NAD  27)  by  NAD  83,  the  readjustment  and 
analysis  activities  involving  European  Datum  1950  (ED  50),  and 
the  availability  of  the  new  Australian  Geodetic  Datum  1984  (AGD 
84)  . 

-  An  extensive  increase  in  the  data  and  types  of  data  needed  to 
develop  an  improved  WGS. 

-  The  availability  of  new  theory  and  techniques  to  support  a  WGS 
improvement  effort. 

WGS  84  has  been  developed  as  a  replacement  for  WGS  72  and 
represents  DMA's  modeling  of  the  earth  from  a  geometric,  geodetic,  and 
gravitational  standpoint  using  data,  techniques,  and  technology 
available  through  early  1984.  It  is  an  improvement  over  WGS  72  in 
several  respects.  New  and  more  extensive  data  sets  and  improved 
computer  software  were  used  in  the  development.  A  more  extensive  file 
of  Doppler-derived  station  coordinates  was  available,  and  for  many  more 
local  geodetic  datums;  improved  sets  of  ground-based  Doppler  and  laser 
satellite  tracking  data  and  surface  gravity  were  available;  and  geoid 
heights  deduced  from  satellite  radar  altimetry  (a  new  data  type)  were 
available  for  oceanic  regions  between  70  degrees  north  and  south 
latitude  (approximately) . 

The  purpose  of  this  publication  i s  to  provide  a  detailed  report  on 
WGS  84  and  its  updates/revisions  which  occurred  since  the  first 
edition.  An  important  feature  of  the  current  edition  includes  the  use 
of  additional  Doppler  and  GPS  survey  information  to  update  the  datum 
transformation  constants. 
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2.  WGS  84  COORDINATE  SYSTEM 


2.1  fifineral 

The  WGS  84  Coordinate  System  is  a  Conventional  Terrestrial 
System  (CTS) ,  realized  by  modifyii.j  the  Navy  Navigation  Satellite 
System  (NNSS) ,  or  TRANSIT,  Doppler  Reference  Frame  (NSWC  9Z-2)  in 
origin  and  scale,  and  rotating  it  to  bring  its  reference  meridian  into 
coincidence  with  the  Bureau  International  de  I'Heure  (BIH) -defined  Zero 
Meridian . 


From  analyses  discussed  in  [1],  it  was  concluded  that  the  NSWC 
9Z-2  Coordinate  System  should  be  modified  by: 

-  Shifting  the  NSWC  9Z-2  origin  by  4 . 5  meters  in  the  negative 
direction  along  the  Z-axis. 

-  Rotating  the  NSWC  9Z-2  Reference  Meridian  (about  the  Z- 
axis)  westward  by  0.814  arc  second  to  the  BIH-defined  Zero 
Meridian  of  1984.0. 

-  Changing  the  NSWC  9Z-2  scale  by  -0.6  x  10“®. 

The  NSWC  9Z-2  Coordinate  System,  modified  in  this  manner  (Table  2.1), 
becomes  (forms)  the  WGS  84  Coordinate  System.  The  origin  and  longitude 
modifications  are  illustrated  in  Figures  2.1  and  2.2,  respectively,  as 
differences  between  NSWC  9Z-2  and  WGS  84.  Use  of  these  modifications 
(Table  2.1)  with  the  Molodensky  Datum  Transformation  Formulas,  after 
modifying  the  formulas  slightly  and  setting  Ax  -  AY  -  0,  provided  the 
A(ji,  AX,,  Ah  formulas  (Table  2.2)  that  produced  the  Doppler  Station  WGS  84 
coordinates  used  to  develop  jLiocal  Geodetic  Datum-to-WGS  84  Datum 
Transformations  (Chapter  7) . 

Thus,  analogous  to  the  BIH-defined  CTS,  or  BIH  Terrestrial 
System  (BTS) ,  the  origin  of  the  WGS  84  Coordinate  System  is  the  center 
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of  mass  of  the  earth;  the  WGS  84  Z-axis  is  in  the  direction  of  the 
Conventional  Terrestrial  Pole  (CTP)  for  polar  motion,  as  defined  by  the 
PIH  for  epoch  1984.0  on  the  basis  of  the  coordinates  adopted  for  the 
BIH  stations;  the  X-axis  is  the  intersection  of  the  WGS  84  reference 
meridian  plane  and  the  plane  of  the  CTP's  equator,  the  reference 
meridian  being  the  Zero  Meridian  defined  by  the  BIH  for  epoch  1984.0  on 
the  basis  of  the  coordinates  adopted  for  the  BIH  stations;  and,  the  Y- 
axis  completes  a  right-handed,  earth-fixed  orthogonal  coordinate  system 
measured  in  the  plane  of  the  above  equator,  90*  east  of  the  X-axis 
(Figure  2.3). 

The  WGS  84  Coordinate  System  origin  and  axes  also  serve  as 
the  geometric  center  and  the  X,  Y,  and  Z  axes  of  the  WGS  84  Ellipsoid. 
(Thus,  the  WGS  84  Coordinate  System  Z-axis  is  the  rotational  axis  of 
the  WGS  84  Ellipsoid.) 

The  WGS  84  Coordinate  System  (reference  frame)  is  the  frame 
of  a  standard  earth  rotating  at  a  constant  rate  around  an  average 
astronomic  pole  (the  CTP) .  However,  the  universe  is  in  motion,  the 
earth  is  nonstandard,  and  events  occur  in  an  instantaneous  world. 
Therefore,  the  WGS  84  Coordinate  System  (CTS)  must  be  related 
mathematically  to  an  Instantaneous  Terrestrial  System  (ITS)  and  to  a 
Conventional  Inertial  System  (CIS) , 

2 . 2  Matnematical  Relationship  Between  the  CIS.  ITS,  and  the 

WGS  84  Coordinate  System 

The  mathematical  relationship  between  the  Conventional 
Inertial  System,  the  Instantaneous  Terrestrial  System,  and  the  WGS  84 
Coordinate  System,  which  is  Identical  to  the  BIH-defined  CTS  in  its 
definition  [1],  can  be  expressed  as; 

CTS  -  [A]  [B]  [C]  [D]  CIS  (2-1) 
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In  Equation  (2-1),  the  rotation  matrices  for  polar  motion 
[A],  sidereal  time  [B] ,  astronomic  nutation  [C] ,  and  precession  [D] 
provide  the  relationship  between  .he  CIS,  defined  by  the  Fundamental 
Katlog  5  (FK5)  System  referenced  to  Epoch  J2000.0  [1],  and  the  WGS  84 
Coordinate  System.  Proceeding  from  right-to-left  in  Equation  (2-1) 
through  matrices  D,  C,  and  B  establishes  the  relationship  between  the 
CIS  and  the  ITS.  Matrix  A  provides  the  relationship  between  the 
Celestial  Ephemeris  Pole  (CEP),  which  approximates  the  instantaneous 
pole  of  the  instantaneous  earth,  and  the  CTP,  or  average  pole  of  the 
standard  earth  associated  with  the  WGS  84  Coordinate  System.  Therefore, 
the  application  of  Matrix  A  completes  the  mathematical  connection 
between  the  WGS  84  Coordinate  System,  an  ECEF  Coordinate  System,  and 
the  CIS,  an  Earth-Centered  Inertial  (ECI)  Cjordinate  System.  For 
detailed  discussion  on  this  subject,  refer  to  [1] . 

Although  tremendous  progress  has  been  made  in  the  last  decade 
in  understanding  and  more  precisely  defining  the  ITS,  the  CTS,  and  the 
CIS,  and  the  mathematical  relationships  between  them  [2],  much  work 
remains  to  be  done.  In  particular,  efforts  to  develop  a  precise 
mathematical  connection  betwvson  stellar  ^optical)  and  radio  (Very  Long 
Baseline  Interferometry,  or  VL)l?I)  systems  and  maintain  the  BlH-defined 
CIS  and  CTS  with  respect  to  a  designated  epoch,  nt^cd  contin-'e. 
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Quantities  Used  in  Converting 
Doppler  System  Coordinates  (NSWC  9Z-2)  to  WGS  84* 
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Navy  Navigatior  satellite  Coordinate  Systera 


WGS  84 


Origin  =  Earth's  center  of  mass 


Z-Axis  =•  The  direction  of  the  Conventional  Terrestrial  Pole  (CTP)  for 
polar  motion,  as  defined  by  the  Bureau  International  de 
I'Heure  (BIH)  on  the  basis  of  the  coordinates  adopted  for 
the  BIH  stations. 

X-Axis  “  Intersection  of  the  WGS  84  Reference  Meridian  Plane  and  the 
plane  of  the  CTP 'a  Equator,  the  Reference  Meridian  being 
the  Zero  Meridian  defined  by  the  BIH  on  the  basis  of  the 
coordinates  adopted  for  the  BIH  stations. 

Y-Axis  =  Completes  a  right-handed.  Earth  Centered,  Earth  Fixed  (ECEF) 

orthogonal  coordinate  system,  measured  in  the  plane  of  the  CTP 
Equator,  90®  East  of  the  X-Axis. 


BIH-Defined  CTP  (1984.0) 


Figure  2.3.  The  WGS  84  Coordinate  System*  Definition 


*  Analogous  to  the  BIH  Defined  Conventional  Terrestrial  System  (CTS) ,  or 
BTS,  1984.0. 
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3.  WGS  84  ELLIPSOID 


3.1  General 

In  geodetic  applications,  three  different  surfaces  or  earth 
figures  are  normally  involved.  In  addition  to  the  earth's  natural  or 
physical  surface,  these  include  a  geometric  or  mathematical  reference 
surface,  the  ellipsoid,  and  an  equipotential  surface  called  the  geoid 
(Chapter  6)  .  In  determining  the  WGS  84  Ellipsoid  and  associated 
parameters,  the  WGS  84  Development  Committee,  in  keeping  with  DMA 
guidance,  decided  quite  early  to  closely  adhere  to  the  thoughts  and 
approach  used  by  the  International  Union  of  Geodesy  and  Geophysics 
(lUGG)  when  the  latter  established  and  adopted  Geodetic  Reference 
System  1980  (GRS  80)  [^].  Accordingly,  a  geocentric  equipotential 

ellipsoid  of  revolution  was  taken  as  the  form  for  the  WGS  84  Ellipsoid. 
The  parameters  selected  to  define  the  WGS  84  Ellipsoid  are  the 

semimajor  axis  (a),  the  earth's  gravitational  constan_^  (GM) ,  the 

normalized  second  degree  zonal  gravitational  coefficient  (  C  20)  the 

angular  velocity  ((0)  of  the  earth  (Table  3.1).  These  parameters  are 

identical  to  those  of  the  GRS  80  Ellipsoid  with  one  minor  exception. 

The  coefficient  form,  used  for  the  second  degree  zonal  is  that  of  the 

WGS  84  Earth  Gravitational  Model  rather  than  the  notation  J2  used  with 

GRS  80.  Accuracy  estimates  (one  sigma)  are  also  included  in  Table  3.1 
for  the  defining  parameters. 

3 . 2  Defining  Parameters 

3.2.1  Semimajor  Axis  (a) 

The  semiraajor  axis  (a)  was  selected  as  one  of  the 
defining  parameters  of  the  WGS  84  Ellipsoid.  Its  adopted  value  and 
estimated  accuracy  (one  sigma)  are: 

a  =  6378137  ±  2  meters.  (3-1) 
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This  value,  which  is  the  same  as  that  of  the  GRS  80  Ellipsoid,  is  two 
meters  (m)  larger  than  the  value  of  6378135  m  adopted  for  the  WGS  72 
Ellipsoid  [4].  As  stated  in  [5],  the  GRS  80,  and  thus  the  WGS  84,  a- 
value  is  based  on  estimates  from  the  1976-1979  time  period,  determined 
using  laser,  Doppler,  radar  altimeter,  laser  plus  radar  altimeter,  and 
Doppler  plus  radar  altimeter  data/techniques.  These  efforts  yielded 
values  from  6378134.5  m  to  6378140  m.  The  best  estimate  was  considered 
to  lie  between  6378135  m  and  6378140  m. 

3.2.2  Earth's  Gravitational  Constant  (GM) 

3. 2. 2.1  GM  With  Earth’s  Atmosphere  Included  (GM). 


The  value  of  the  earth's  gravitational  constant, 
adopted  as  one  of  the  four  defining  parameters  of  the  WGS  84  Ellipsoid, 
and  its  one-sigma  accuracy  estimate  are: 

GM  -  (3986005  ±0.6)  x  lO^m^s’^  .  (3-2) 


This  value  includes  the  mass  of  the  atmosphere  and  is  based  on  several 
types  of  space  measurements.  These  measurement  types  and  the  associated 
estimates  for  GM  are  [3]  : 


Spacecraft  radio  tracking  .  (3986005.0  ±0.5) 

Lunar  laser  data  analysis  .  (3986004.6  ±0.3) 

Satellite  laser  range  measurements  .  (3986004.4  ±0.2) 


From  these  results,  the  representative  value  in  Equation 
consistent  with  the  data  used,  was  then  adopted. 


X  10®m^s"^ 

X  10®m^s"^ 

X  10“m^s"^ 

(3-2)  for  GM, 
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3. 2. 2. 2  GM  of  the  Earth’s  Atmosphere 


For  some  applications,  it  is  necessary  to  either 
have  a  GM  value  for  the  earth  which  does  not  include  the  mass  of  the 
earth's  atmosphere,  or  have  a  GM  value  for  the  earth's  atmosphere 
itself.  For  this,  it  is  necessary  to  know  both  the  mass  of  the  earth's 
atmosphere,  M;^,  and  the  universal  gravitational  constant,  G. 

Using  the  value  recommended  for  G  [6]  by  the 
International  Association  of  Geodesy  (lAG)  ,  and  the  more  recent  value 
for  M;^  [7],  the  product  GM^^  to  two  significant  digits  yields  the  value 

currently  recommended  by  the  lAG  for  this  constant  [6]  .  This  value, 
with  an  assigned  accuracy  estimate,  was  adopted  for  use  with  WGS  84 : 

GM;^  =  (3.5  ±  0.1)  X  10®m'^s-2  (3-3) 

3. 2. 2. 3  GM  With  J:arth's  Atmosphere  Excluded  (GM'I 

The  earth's  gravitational  constant  with  the  mass 
of  the  earth's  atmosphere  excluded  (GM'),  was  obtained  by  subtracting 
GM;^,  Equation  (3-3),  from  GM,  Equation  (3-2) 

GM'  =  (3986001.5  ±  0.6)  x  lO^m^s"^  (3-4) 

The  fact  that  the  WGS  84  value  for  GM',  Equation  (3-4),  is  given  to  one 
more  digit  than  the  WGS  84  value  for  GM,  Equation  (3-2),  does  not  imply 
that  GM'  is  known  more  accurately  than  GM.  The  additional  digit  used 
with  GM'  only  reflects  a  desire  to  maintain  consistency  between  the 

various  WGS  84  parameters  and  correction  terms.  In  fact,  GM*  is  known 
less  well,  due  to  the  uncertainty  introduced  via  GM;^ .  The  lack  of  a 

more  realistic  accuracy  value  for  GM^  prevents  acknowledgment  of  this 
in  the  above  one-sigma  accuracy  estimate  for  GM' . 


3-3 


3.2.3  Normalized  S^ond  Degree  Zonal  Gravitational 
Coefficient  C  20 

Another  defining  parameter  of  the  WGS  84  Ellipsoid  is 
the  normalized  second  degree  zonal  gravitational  coefficient,  C  20^ 

which  has  the  following  value  and  assigned  accuracy  (one  sigma) : 

C  20  =  (-484.16685  +  0.00130)  x  10"^  (3-5) 

This  C  20  value  was  obtained  from  the  adopted  GRS  80  value  for  J2  [3], 

(  J2°“  J20  )  ' 


J2  =  108263  X  10'® 


(3-6) 


by  using  the  mathematical  relationship 

"^20  =  -J2/(5)1''2  (3-7) 

and  truncating  the  result  to  eight  significant  digits. 

In  keeping  with  the  GRS  80  value  for  J2,  the  C  20  value 

for  the  WGS  84  Ellipsoid  also  does  not  include  the  permanent  tidal 
deformation.  This  effect,  usually  represented  by  6J2,  is  due  to  the 

attraction  of  the  earth  by  the  sun  and  moon.  It  has  the  magnitude  [8]: 

5j.  =  9.3  X  10'®  (3-8) 


or,  equivalently 


6  C  20  =  -4.16  X  10-5.  (3_9) 

This  quantity  would  be  added  to  C  20/  Equation  (3-5)  ,  if  it  were 
desired  to  have  C 20  include  the  permanent  tidal  deformation. 
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3.2.4  Aniaular  Velocity  of  the  Earth_(a)) 


The  value  of  (O  used  as  one  of  the  defining  parameters 
of  the  WGS  84  (and  GRS  80)  Ellipsoid  and  its  accuracy  estimate  (one 
sigma)  are: 

CD  =  (7292115  ±  0.1500)  x  10'“  radians/second  (S-IO; 

This  value,  for  a  standard  earth  rotating  with  a  constant  angular 
velocity,  is  an  lAG  adopted  value  for  the  true  angular  velocity  of  the 
earth  which  fluctuates  with  time.  However,  for  most  geodetic 
applications  which  require  angular  velocity,  these  fluctuations  do  not 
have  to  be  considered. 

Although  CO  is  suitable  for  use  w'ith  a  standard  earth 

and  the  WGS  84  Ellipsoid,  it  is  the  International  Astronomical  Union 
(lAU),  or  the  GRS  67,  version  of  this  value  (CD') 

CD'  =  7292115.1467  x  10"^^  radians/second  (3-11) 

that  was  used  with  the  new  definition  of  time  [9] . 

For  consistent  satellite  applications,  the  value  of  the 
earth's  angular  velocity  (O)')  from  Equation  (3-11),  rather  than  CD, 

should  be  used  in  the  formula 

CD*  =  CD'+  m  (3-12) 

to  obtain  the  angular  velocity  of  the  earth  in  a  precessing  reference 
frame  (CD*)  .  In  the  above  equation  [9]  [10]  : 

m  =  precession  rate  in  right  ascension 

m  =  (7.086  X  10"^^  +  4.3  x  10'^^  T^)  radians/second 

(3-13) 


3-5 


Tu  “  Julian  Centuries  from  Epoch  J2000.0 

Tu  ■■  du/36525 

du  “  Nuniber  of  days  of  Universal  Time  (UT)  from  Julian 

Date  (JD)  2451545.0  UTl,  taking  on  values  of  +  0.5, 

i  1.5,  i  2.5,... 

du  -  JD  -  2451545. 

Therefore,  the  angular  velocity  of  the  earth  in  a 
precessing  reference  frame,  for  satellite  applications,  is  given  by: 

(0*  -  (7292115.8553  x  lO'^  +  4.3  x  T,;) 

radians/second  (3-14) 

3.3  Derived  Geometric  and  Physical  Constants 

3.3.1  General 

Many  parameters  associated  with  the  WGS  84  Ellipsoid, 
other  than  the  four  defining  parameters  (Table  3.1),  are  needed  for 
geodetic  and  gravimetric  applications.  Using  the  four  defining 
parameters,  it  is  possible  to  derive  these  associated  constants.  The 
more  commonly  used  geometric  and  physical  constants  associated  with  the 
WGS  84  Ellipsoid  are  listed  in  Tables  3.2  and  3.3.  The  formulas  used  in 
the  calculation  of  these  constants  are  primarily  from  [3]  and  [11] . 

The  defining  parameters  are  considered  to  be  exact .  On 
the  other  hand,  the  other  constants  are  derived.  Users  are  reminded 
that  the  derived  constants  must  retain  the  listed  significant  digits  if 
consistency  between  the  magnitudes  of  the  various  parameters  is  to  be 
maintained.  These  constants  should  always  be  calculated  to,  and  used 
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with/  the  number  of  digits  required  to  maintain  the  consistency  needed 
for  each  specific  application. 

3.3.2  Relevant  Miscellaneous  Constants/Conversion  Factors 

In  addition  to  the  four  defining  parameters  of  the  WGS 
84  Ellipsoid  (Table  3.1),  necessary  for  describing  (representing)  the 
ellipsoid  geometrically  and  gravimetrically,  and  the  derived  sets  of 
commonly  used  geometric  and  physical  constants  associated  with  the  WGS 
84  Ellipsoid  (Tables  3.2  and  3.3),  two  other  important  constants  are  an 
integral  part  of  the  definition  of  WGS  84.  These  constants  are  the 
velocity  of  light  (c)  and  the  dynamical  ellipticity  (H) . 

The  currently  accepted  value  for  the  velocity  of  light 
in  a  vacuum  (c)  is  [12]  : 

c  -  (299792458  ±  1.2)  m  s"^  .  (3-15) 

This  value  is  officially  recognized  by  both  the  lAG  [6]  and  lAU  [10], 
and  has  been  adopted  for  use  with  WGS  84 . 

The  dynamical  ellipticity  (H)  is  necessary  for 
determining  the  earth's  principal  moments  of  inertia,  A,  B,  and  C.  In 
the  literature,  H  is  variously  referred  to  as  dynamical  ellipticity, 
mechanical  ellipticity,  or  the  processional  constant .  It  is  a  factor  in 
the  theoretical  value  of  the  rate  of  precession  of  the  equinoxes,  which 
is  well  known  from  observation.  In  a  recent  lAG  report  on  fundamental 
geodetic  constants  [8],  the  following  value  for  the  reciprocal  of  H  was 
given  in  the  discussion  of  moments  of  inertia: 

1/H  -  305.4413  ±  0.0005.  (3-16) 

For  consistency,  this  value  has  been  adopted  for  use  with  WGS  84 . 
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Values  of  the  velocity  of  light  in  a  vacuum  and  the 
dynamical  ellipticity  adopted  for  use  with  WGS  84  are  listed  in  Table 
3.4  along  with  other  WGS  84  associated  constants  used  in  special 
applications.  Factors  for  effecting  a  conversion  between  meters,  feet, 
and/or  nautical  and  statute  miles  are  also  given  in  the  table. 

3.4  Comments 

The  four  defining  parameters  <  a,  C  20f  to,  GM  )  of  the  WGS  84 
Ellipsoid  were  used  to  calculate  the  more  commonly  used  geometric  and 
physical  constants  associated  with  the  WGS  84  Ellipsoid.  As  a  result  of 
the  use  of  C  20  if^  the  form  described,  the  derived  WGS  84  Ellipsoid 

parameters  are  slightly  different  from  their  GRS  80  Ellipsoid 
counterparts.  Although  these  minute  parameter  differences  and  the 
conversion  of  the  GRS  80  Jj-value  to  C  20  ^^6  insignificant  from  a 

practical  standpoint,  it  has  been  more  appropriate  to  refer  to  the 
ellipsoid  used  with  WGS  84  as  the  WGS  84  Ellipsoid. 

_  In  contrast,  since  NAD  83  does  not  have  an  associated  EGM,  the 

J2  to  C  20  conversion  does  not  arise  and  the  ellipsoid  used  with  NAD  83 

by  the  National  Geodetic  Survey  (NGS)  is,  in  name  and  in  both  defined 
and  derived  parameters,  the  GRS  80  Ellipsoid.  Although  it  is  important 
to  know  that  these  small  undesirable  inconsistencies  exist  between  the 
WGS  84  and  GRS  80  Ellipsoids,  from  a  practical  application  standpoint 
they  are  insignificant.  This  is  especially  true  with  respect  to  the 
defining  parameters.  Therefore,  as  long  as  the  preceding  is  recognized, 
it  can  be  stated  that  WGS  84  and  NAD  83  are  based  on  the  same 
ellipsoid. 
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Parameters  1  Notation  I  Magnitude  Accuracy  (ICT) 
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Table  3.4 

Relevant  Miscellaneous  Constants 
and  Conversion  Factors 


Constant 

Symbol 

1 

Numerical  Value 

Velocity  of  Light 

c 

299792458  m  s'^ 

(in  a  Vacuum) 

Dynamical  Ellipticity 

H 

1/305.4413 

Earth's  Angular  Velocity  [for 
Satellite  Applications;  see 

CD* 

(7292115.8553  x  10'“ 

+  4.3  X  10'^^  Tu)rad  s'^ 

Equation  (3-14)] 

Universal  Constant  of  Gravitation 

G 

6.673  X  10"“m^s'^)cg'^ 

GM  of  the  Earth's  Atmosphere 

GM^ 

3.5  X  10®m^s'^ 

Earth's  Gravitational  Constant 

GM' 

3986001.5  X  10®m^s'2 

(Excluding  the  Mass  of  the 

Earth's  Atmosphere) 

Earth's  Principal 

A 

8.0091029  X  10^^  leg  m^ 

Moments  of  Inertia 

B 

8.0092559  x  10^''  kg  m^ 

(Dynamic  Solution) 

C 

8.0354872  x  lO^’^  kg  m^ 

Conversion 

Factors 

1  Meter 
1  Meter 

1  International  Foot 
1  US  Survey  Foot 
1  US  Survey  Foot 


=  3.28083333333  US  Survey  Feet 
=  3.28083989501  International  Feet 
=  0.3048  Meter  (Exact) 

=  1200/3937  Meter  (Exact) 

=  0.30480060960  Meter 


=  1852  Meters  (Exact) 

1  International  Nautical  Mile  =  6076.10333333  US  Survey  Feet 

=  6076.11548556  International  Feet 


1  International  Statute  Mile  =  1609.344  Meters  (Exact) 

=  5280  International  Feet  (Exact) 


Tu  “  Julian  Centuries  from  Epoch  J2000.0 


4.  WGS  84  ELLIPSOIDAL  GRAVITY  FORMULA 


4.1  General 

In  Section  3.1,  the  WGS  84  Ellipsoid  is  identifitid  as  being  a 
geocentric  equipotential  ellipsoid  of  revolution.  An  equipotential 
ellipsoid  is  simply  an  ellipsoid  defined  to  be  an  equipotential 
surface,  i.e.,  a  surface  on  which  all  values  of  the  gravity  potential 
are  equal.  Given  an  ellipsoid  of  revolution,  it  can  be  made  an 
equipotential  surface  of  a  certain  potential  function,  the  theoretical 
(normal)  gravity  potential  (U) .  This  theoretical  gravity  potential  can 
be  uniquely  determined,  independent  of  the  density  distribution  within 
the  ellipsoid,  by  using  any  system  of  four  independent  constants  as  the 
defining  parameters  of  the  ellipsoid.  As  noted  earlier  for  the  WGS  84 
Ellipsoid  (Chapter  3),  these  are  the  semimajor  axis  Ja) ,  the  normalized 
second  degree  zonal  gravitational  coefficient  (  C 20) <  the  earth's 

angular  velocity  ((O)  ,  and  the  earth's  gravitational  constant  (GM)  . 

Theoretical  gravity  (y) ,  the  gradient  of  U,  is  given  on  (at) 
the  surface  of  the  ellipsoid  by  the  closed  formula  of  Somigliana  [13] : 

Y  “  Yo  cos^<t)  +  b  Yp  sin^<|>) )  /  (a^cos^<|)  +  b^sin^<j))  (4'-l) 


where 


a,  b  "  semimajor  and  semiminor  axes  of  the  ellipsoid, 
respectively 

Ye'Yp  "  theoretical  gravity  at  the  equator  and  poles, 
respectively 

<))  •»  geodetic  latitude. 
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Thus,  the  equipotential  ellipsoid  serves  not  only  as  the  reference 
surface  or  geometric  figure  of  the  earth,  but  leads  to  a  closed  formula 
for  theoretical  gravity  at  the  ellipsoidal  surface. 

4.2  Analytical  and  Numerical  Forms 

The  closed  gravity  formula  of  Somigliana  in  the  form  [3] 

7  “  Ye  (i  +  k  sin^<t>)/{l-  e^sin^(|)) (4-2) 

has  been  selected  as  the  official  VJGS  84  Ellipsoidal  Gravity  Formula. 

In  Equation  (4-2)  ; 

k  =  (b  Yp/a  Ye)  “  ^  (4-3) 

e^  «=  square  of  the  first  eccentricity  of  the  ellipsoid. 

Equation  (4-2)  was  selected  for  use  with  WGS  84  in  preference  to 

Equation  (4-1)  since  it  is  more  convenient  for  numerical  computations 
and  explicitly  contains  only  Ye  the  first  factor  in  the  equation. 

The  analytical  and  numerical  forms  of  the  WGS  84  Ellipsoidal 
Gravity  Formula  are  provided  in  Table  4.1. 
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Table  4 . 1 


WGS  84 

Ellipsoidal  Gravity  Formula 

Provides  Gravity  Values  at  (on)  the  Surface  of  the  WGS  84  Ellipsoid 

Notation 

Y  >»  Acceleration  of  a  unit  test  mass  due  to  theoretical  gravity. 

Yo  "  Acceleration  at  the  equator  (on  the  WGS  84  Ellipsoid)  of  a  unit 
test  mass  due  to  theoretical  gravity. 

k  -  Consr.ant  =  (b  Yp/^  Ye)  “  1 

a  “  Semimajor  axis  (WGS  84  Ellipsoid) 
b  =  Semiminor  axis  (WGS  84  Ellipsoid) 

Yp  =  Theoretical  gravity  at  the  poles  (on  the  WGS  84  Ellipsoid) 
(j)  «=  Geodetic  latitude 

e^  =  First  eccentricity  squared  (WGS  84  Ellipsoid) . 

Malytical  Tocm 

Y  =  Ye  (1  +  k  sin^^)/(l-  e^sin^(j)) 

Numerical  Form 

Y  «=  978032.67714  (1  +  0.00193185138639  sin^}  / 

(1  -  0.00669437999013  sin^(|>)^^^  x  10”^  m/second^  (mgal) 

An  acceleration  due  to  gravity  of  1  x  10“^  m/second^  =  1  mgal 
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5.  WGS  84  GRAVITY  MODELING 


5.1  Earth  Gravitational  Model _ (EGM^ 

The  form  of  the  WGS  84  EGM  is  a  spherical  harmonic  expansion 
(Table  5.1)  of  the  gravitational  potential  (V).  The  WGS  84  EGM, 
complete  through  degree  (n)  and  order  (m>  180,  is  comprised  of  32755 
coefficients . 

The  coefficients  through  n=m=41  were  obtained  from  a  weighted 
least  squares  solution  of  a  normal  equation  matrix  developed  by 
combining  individual  normal  equation  matrices  formed  from  Doppler 
satellite  tracking  data,  satellite  laser  ranging  data,  surface  gravity 
data,  oceanic  geoid  heights  deduced  from  satellite  radar  altimeter 
data,  Navstar  Global  Positioning  System  (GPS)  data,  and  "lumped 
coefficients".  The  effect  (contribution)  of  coefficients  through  n="m=41 
was  removed  from  a  worldwide  l°xl°  mean  gravity  anomaly  field  leaving  a 
worldwide  residual  l°xl°  mean  gravity  anomaly  field.  The  WGS  84  EGM 
coefficients  from  n=42,  m=0  through  n=m=180  were  then  determined 
independently  via  harmonic  analysis  using  the  residual  field.  The 
coefficients  through  n=m=41  from  the  weighted  least  squares  solution 
and  the  coefficients  above  n=m=41  from  the  independent  harmonic 
analysis  comprise  the  n=m=180  WGS  84  EGM. 

The  WGS  84  EGM  through  n=m=180  is  to  be  used  when  calculating 
WGS  84  Geoid  Heights,  WGS  84  gravity  disturbance  components  (or 
deflection  of  the  vertical  components),  and  WGS  84  l°xl°  mean  gravity 
anomalies  via  spherical  harmonic  expansions.  Expansions  to  this  degree 
and  order  (n=m=180)  are  needed  to  accurately  model  variations  in  the 
earth's  gravitational  field  on  or  near  the  earth's  surface. 

The  WGS  84  EGM  through  n=m=41  is  more  appropriate  for 
satellite  orbit  calculation  and  prediction  purposes.  The  use  of  higher 
degree  and  order  models  for  such  applications  is  not  recommended  at 
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this  time.  However,  if  required  for  a  special  application,  DMA  and 
other  DoD  users  will  need  to  conduct  orbital  analyses  and  ascertain  the 
EGM  truncation  level  that  is  "best”  suited  for  the  satellite  project 
involved. 


The  WGS  84  EGM  through  n=m=180  is  available  on  magnetic  tape 
in  normalized  form.  The  WGS  84  EGM  through  n=m=41  is  available  on  a 
separate  magnetic  tape  in  both  normalized  and  conventional  form. 
However,  the  WGS  84  EGM  coefficients  through  n=m=18  are  provided  in 
Table  5.2  in  normalized  form. 

Accuracy  values  are  not  available  for  all  of  the  WGS  84  EGM 
coefficients.  However,  an  error  covariance  matrix  is  available  for 
those  coefficients  through  n=m=41  determined  from  the  weighted  least 
squares  solution.  Gravity  anomaly  degree  variances  are  given  in  Table 
5.3  for  the  WGS  84  EGM  (n=m=180) .  Requesters  having  a  need  for  the  full 
WGS  84  EGM  and/or  its  error  data  should  forward  their  correspondence  to 
the  address  listed  in  the  PREFACE. 

5.2  Gravity  Potential  (W) 

Using  the  WGS  84  EGM  model  (Table  5.1),  the  earth's  total 
gravity  potential  (W)  is  then  defined  as 


W  =  V  +  <[)  (5-1) 

where  O  is  the  centrifugal  potential  due  to  the  earth’s  rotation.  If 
CO  is  the  angular  velocity  [Equation  (3-10)],  then 

O  =  I  (02(x2  +  Y^)  (5-2) 

where  X  and  Y  are  the  geocentric  coordinates  of  the  rotating  mass  in 
the  WGS  84  reference  frame  (See  Figure  2.3). 
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Table  5.1 


Form  of  the  WGS  84 
Earth  Gravitational  Model 


V  =  ^  [  1  +  Z" 

r  L  -L  rr 

-0  "  Pnm(sin(})')  (Cn„cosm>.+  Sn^sinmX)] 

Parameter 

Definition 

V 

=  Gravitational  potential  function  (m^s'^) 

GM 

“  Earth's  gravitational  constant 

r 

”  Radius  vector  from  the  earth's  center  of 

mass 

a 

“  Semiraajor  axis  of  the  WGS  84  Ellipsoid 

n,m 

-  Degree  and  order,  respectively 

r 

=»  Geocentric  latitude** 

X 

=•  Geocentric  longitude  =  geodetic  longitude** 

^  nm'  ^  nm 

-  Normalized  gravitational  coefficients* 

P  nm(sin  (j)'  ) 

“  Normalized  associated  Legendre  function 

..  .  .  _  1/2 

f  (n-m)  !  {2n+l)k1  , 

-  L-' - '(n+m)!  J  f  > 

£’nm(sin  4)') 

=  Associated  Legendre  function 

*  See  next  page . 

**  Latitude  is  positive  north  and  longitude  is  positive  east  (0°  to 
180°) 
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Table  5.1  (Cont'd) 


Form  of  the  WGS  84 
Earth  Gravitational  Model 


Parameter 

Definition 

P„„(sin  (t)') 

=  (cos  (j)')"’  ,  [Pn(sin  ()>')] 

d (sin  (])')"' 

Pn  (sin  (()' ) 

=  Legendre  polynomial 

Pn  (sin  (()' ) 

*=  (sin^(j)'  -  1)" 

2"n!  d(sin  (j)')" 

*Note: 


(n+m)  ! 


1/2 


(n-m) !  (2n+l)  k 


where 


=  Conventional  gravitational  coefficients 

For  m=0,  k=l/ 
m>l,  k=2 . 
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Table  5 . 3 


WGS  84  EGM 

Gravity  Anomaly  Degree  Variances  (c^) 


Degree 

Degree 

Variances 

Degree 

Degree 

Variances 

Degree 

Degree 

Variances 

2 

7.6 

41 

2.8 

80 

2.3 

3 

33.9 

42 

2.7 

81 

2.6 

4 

19.2 

43 

2.4 

82 

2.8 

5 

21.0 

44 

2.8 

83 

2.7 

6 

19.4 

45 

2.7 

84 

2.4 

7 

19.3 

46 

2.7 

85 

2.2 

8 

10.9 

47 

3.1 

86 

2.7 

9 

11.5 

48 

2.5 

87 

2.4 

10 

9.7 

49 

2.5 

88 

2.2 

11 

6.4 

50 

2.8 

89 

2.0 

12 

2.6 

51 

2.8 

90 

2.0 

13 

7.4 

52 

2.6 

91 

2.3 

14 

3.2 

53 

3.1 

92 

2.0 

15 

3.4 

54 

2.8 

93 

2.2 

16 

3.9 

55 

3.0 

94 

2.0 

17 

3.6 

56 

3.0 

95 

1.8 

18 

3.6 

57 

3.0 

96 

2.0 

19 

3.3 

58 

2.6 

97 

1.9 

20 

3.1 

59 

3.0 

98 

2.1 

21 

3.2 

60 

2.6 

99 

1.7 

22 

3.6 

61 

2.5 

100 

1.8 

23 

2.7 

62 

2.9 

101 

1.7 

24 

2.6 

63 

2.5 

102 

2.1 

25 

2.9 

64 

2.7 

103 

2.3 

26 

2.4 

65 

2.1 

104 

1.8 

27 

1.9 

66 

2.6 

105 

1.8 

28 

2.4 

67 

2.5 

106 

1.7 

29 

2.4 

68 

2.6 

107 

1.8 

30 

2.8 

69 

2.9 

108 

1.9 

31 

2.9 

70 

2.3 

109 

2.0 

32 

4.1 

71 

2.3 

110 

2.0 

33 

3.4 

72 

2.7 

111 

1.7 

34 

5.0 

73 

2.4 

112 

1.6 

35 

4.4 

74 

2.6 

113 

1.8 

36 

3.6 

75 

2.2 

114 

1 . 6 

37 

3.4 

76 

2.3 

115 

1.9 

38 

2.8 

77 

2.3 

116 

1.7 

39 

3.5 

78 

2.4 

117 

1.7 

40 

3.6 

79 

2.2 

118 

1.6 

Units  =  (1  X  10"^  m/second^) ^  or  mgal^ 


*  See  next  page . 
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Table  5.3  (Cont'd) 


WGS  84  EGM 

Gravity  Anomaly  Degree  Variances  (c„)  * 


Degree 

Degree 

Variances 

Degree 

Degree 

Variances 

Degree 

Degree 

Variances 

mM 

1.6 

140 

1.3 

161 

0.9 

1.6 

141 

1.1 

162 

0.8 

1.5 

142 

1.0 

163 

0.9 

122 

1.3 

143 

1.1 

164 

1.0 

123 

1.6 

144 

1.1 

165 

0.9 

124 

1.6 

145 

1.0 

166 

0.8 

125 

1.5 

146 

1.0 

167 

0.9 

126 

1.3 

147 

1.0 

168 

0.8 

127 

1.5 

148 

1.2 

169 

0.8 

128 

1.2 

149 

1.0 

170 

0.8 

129 

1.4 

150 

1.0 

171 

0.9 

130 

1.3 

151 

1.0 

172 

0.7 

131 

1.3 

152 

1.0 

173 

0.8 

132 

1.4 

153 

1.0 

174 

0.8 

133 

1.4 

154 

1.0 

175 

0.8 

134 

1.2 

155 

0.9 

176 

0.8 

135 

1.2 

156 

1.0 

177 

0.6 

136 

1.2 

157 

0.8 

178 

0.7 

137 

1.3 

158 

0.8 

179 

0.8 

138 

1.3 

159 

0.8 

180 

0.8 

139 

1.2 

160 

0.9 

Units  =  (1  X  10"^  m/second^)^  or  mgal^ 


*Formula  for  computing  gravity  anomaly  degree  variances  (c^) 

c„  -  (  c'",  +  'sIj 

Cn  =  Gravity  anomaly  degree  variance  in  mgal^ 

_  for  degree  n 

Y  “  Average  value  of  theoretical  gravity 

Y  =  979764.46561  mgal  (based  on  WGS  84 

Ellipsoidal  Gravity  Formula) 

C  =  Normalized  gravitational  coefficients  of 
S  nn,  degree  n  and  order  m 
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6.  WGS  84  GEO ID 


6.1  General 

In  geodetic  applications,  three  different  surfaces  or  earth 
figures  are  normally  involved.  In  addition  to  the  earth's  natural  or 
actual  surface,  the  other  two  include  a  geometric  or  mathematical 
figure  taken  to  be  an  equipotential  ellipsoid  of  revolution  (Chapter 
3)  ,  and  a  physical  figure  defined  as  an  equipotential  surface  in  the 
earth's  gravity  field,  or  the  geoid. 

The  equation 

W(X,Y,Z)  =  constant  (6-1) 

defines  a  family  of  equipotential  surfaces  (GEOPS)  of  the  earth's 
gravity  field.  The  geoid  is  that  particular  equipotential  GEOP  for 
which  the  constant  in  Equation  (6-1)  is  equal  to  Wq  (or  the  ellipsoidal 
potential  Uq)  defined  in  Table  3.3. 

For  some  practical  applications,  the  geoid,  defined  as  above, 
is  approximated  by  the  mean  sea  level  (msl)  over  the  oceans  (or  its 
hypothetical  extension  under  the  land  masses) .  It  may  be  necessary  to 
clarify  that  the  msl  is  not  an  equipotential  surface  and  in  its 
simplest  definition  would  comprise  a  mean  of  sea  level  surfaces 
approximated  and  observed  over  18.67  years. 

In  a  mathem.atical  sense,  the  geoid  is  also  defined  (or 
realized)  as  so  many  meters  above  (+N)  or  below  (-N)  the  ellipsoid,  the 
geometric  figure. 

In  the  definition  of  the  geoid,  a  great  practical  importance 
exists:  the  geoid  can  serve  as  the  approximation  for  the  vertical  datum 
reference  surface  for  mean  sea  level  heights  (H)  .  In  areas  where 

*  Note  the  other  usage  of  this  symbol  in  Equation  (3-16)  and  Table  3.4. 
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general  elevation  data  is  not  available  from  conventional  leveling,  the 
"approximate"  determination  of  the  H-values  can  be  achieved  using  the 
equation 

h  ~  N  +  H  (6-2) 

H  ~  h  -  N  (6-3) 

where : 

h  =  geodetic  height  =  height  relative  to  the  ellipsoid 
N  =  geoid  height 

H  =  orthometric  height  relative  to  the  geoid  (or, 
approximately,  mean  sea  level  height) 


Equation  (6-3)  illustrates  the  use  of  geoid  heights  in  the 
determination  of  H-values  from  geodetic  heights  derived  using  satellite 
positions  (e.g.,  TRANSIT  or  Navstar  Global  Positioning  System)  located 
on  the  earth's  physical  surface  or  aboard  a  vehicle  operating  near  the 
earth's  surface. 


6 . 2  Formulas  and  Representations/Analvsis. 


6.2.1  Formulas 


The  WGS  84  Geoid  Heights  were  calculated  using  a 
spherical  harmonic  expansion  and  the  WGS  84  EGM  coefficients  through 
n=m=180.  The  formula  for  calculating  WGS  84  Geoid  Heights  has  the  form: 


N  = 


GM 

ry 


r  y'^max 

L^n=2 


‘^m=0 


(i) 


(  c  ^„cos 


mX  + 


S  nmSin 


mX)  P  nm(sin  (j)'  ) 


(6-4) 


v/here  N  is  the  geoid  height  in  meters,  y  is  theoretical  gravity 
calculated  using  the  WGS  84  Ellipsoidal  Gravity  Formula  (Table  4.1), 
and  all  other  quantities  are  defined  as  for  the  WGS  84  EGM  with  one 
exception.  In  Equation  (6-4),  the  even  degree  zonal  coefficients  of 
subscripts  2  through  10  are  coefficient  differences  between  the  WGS  84 
EGM  minus  normalized  gravity.  (See  Sections  6.2.1  and  6.2.2  in  [1]) . 
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6.2.2  Representations/Analysia 


The  geoid  can  be  depicted  as  a  contour  chart  which 
shows  the  deviations  of  the  geoid  from  the  ellipsoid  selected  as  the 
mathematical  figure  of  the  earth.  Figure  6.1  is  a  worldwide  WGS  84 
Geoid  Height  Contour  Chart  developed  from  a  worldwide  l°xl°  grid  of 
geoid  heights  calculated  by  using  WGS  84  numerical  data  and  the  WGS  84 
EGM  coefficients  through  n=m=18  in  Equation  (6-4) . 

Table  6.1  contains  a  worldwide  10°xl0°  grid  of  WGS  84 
Geoid  Heights  calculated  using  WGS  84  EGM  coefficients  through  n=m“180. 

A  worldwide  l°xl°  grid  of  WGS  84  Geoid  Heights  was 
computed  and  compared  with  a  similar  grid  of  WGS  72  Geoid  Heights 
referenced  to  the  WGS  72  Ellipsoid.  The  root-mean-square  (RMS) 
difference  was  +4.6  meters,  with  the  largest  positive  and  negative 
differences  being  24  and  -23.5  meters,  respectively.  Of  the  64,800 
geoid  height  differences,  21.36  percent  (13,841  differences)  were 
larger  then  5  meters . 

The  RMS  WGS  84  Geoid  Height,  taken  worldwide  on  the 
basis  of  a  l°xl°  grid,  is  30.5  meters.  This  RMS  value  indicates  how 
well  the  WGS  84  Ellipsoid,  taken  as  the  mathematical  figure  of  the 
earth,  fits  the  earth's  geoidal  surface. 

The  WGS  84  Geoid  Heights  have  an  error  range  of  +2  to 
±6  meters  (one  sigma) ,  and  are  known  to  accuracies  of  ±2  to  ±3  meters 
over  approximately  55  percent  of  the  earth.  Approximately  93  percent  of 
the  earth  has  WGS  84  Geoid  Heights  of  accuracy  better  th^n  +4  meters. 

6.3  Availability  of  WGS  84  Geoid  Height  Data 

The  WGS  84  Geoid  Heights,  the  related  data,  and  products, 
which  can  be  provided  to  requesters  (see  PREFACE),  are: 

-  A  worldwide  WGS  84  Geoid  Height  Contour  Chart  with  5  meters 
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contour  interval.  If  needed,  contour  charts  of  other  physical  sizes, 
geographic  areas,  contour  intervals  and  scales  can  be  provided. 

-  A  magnetic  tape  containing  the  worldwide  l°xl“  or  30'x30' 
grid  of  WGS  84  Geoid  Heights. 

-  A  Bi-Linear  Interpolation  program  (Table  6.2)  for 
interpolation  of  WGS  84  Geoid  Heights  at  random  points.  Users  are 
advised  to  check  interpolation  error(s)  pertaining  to  application (s) . 

-  A  Computer  Program  for  computation  of  WGS  84  Geoid  Heights 
at  a  specified  grid  interval  or  at  random  points  with  associated 
documentation  and  appropriate  test  cases. 
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Table  6.2 


Interpolation  of  WGS  84  Geoid  Heights 
Bi-Linear  Interpolation  Method  (Formula) 

Np((|)A)  =  ao  +  aiX  +  ajY  +  a3XY 

-  Geoid  height  (N)  to  be  interpolated  at 
Point  P  (())/  X) 

Sq  “  N3 

ai  “  Nj  -  N3 

aj  -  N4  -  Ni 

aa  -  Ni  +  N3  -  N2  -  N4 

X  =  (X  -  li)/[\2  -  K) 

Y  -  -  ^i)  /  (<j>2  -  <{>i) 

(j)  -  Geodetic  latitude  of  Point  P 

X  -  Geodetic  longitude  of  Point  P 

Nj^,  Nj,  N3,  N4  =“  Known  geoid  heights  at  grid  points  used  in  the 
interpolation  process 

[See  Associated  Coordinate  System,  Figure  6.2] 

Note:  Use  of  consistent  units  is  necessary. 
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7.  WGS  84  RELATIONSHIPS  WITH  OTHER  GEODETIC  SYSTEMS 


7.1  General 

One  of  the  principal  purposes  of  a  world  geodetic  system  is  to 
provxde  the  means  whereby  local  geodetic  datums  can  be  referenced  to  a 
single  geocentric  system.  The  number  of  local  geodetic  datums,  or  local 
horizontal  datums,  requiring  such  referencing  is  extensive.  Counting 
island  and/or  astronomic-based  datums,  the  number  exceeds  several 
hundred.  To  accomplish  the  conversion,  local  geodetic  datum  and  WGS 
coordinates  are  both  required  at  one  or  more  sites  within  the  local 
datum  area  so  that  local  geodetic  datum-to-WGS  datum  shifts  can  be 
formed.  Satellite  stations  positioned  within  WGS  84,  and  with  known 
local  geodetic  datum  coordinates,  were  the  basic  ingredients  in  the 
development  of  local  geodetic  datum-to-WGS  84  datum  shifts. 

The  most  accurate  approach  for  obtaining  WGS  84  coordinates  is 
to  acquire  satellite  tracking  data  at  the  site  of  interest  and  position 
it  directly  in  WGS  84  using  the  Satellite  Point  Positioning  technique 
[1] .  However,  it  is  unrealistic  to  presume  that  use  of  this  technique 
will  always  be  possible.  In  such  cases,  the  transformation  from  WGS  72 
to  WGS  84  or  from  local  geodetic  daturas  to  WGS  84  should  be  used 
(sections  7.2  and  7.3). 

7.2  WGS  72-to-WGS  84  Transformation 

Situations  arise  where  only  WGS  72  coordinates  are  available 
for  a  site.  In  such  instances,  the  WGS  72-to-WGS  84  Transformation 
listed  in  Table  7.1  can  be  used  with  the  following  equations  to  obtain 
WGS  84  coordinates  for  the  sites: 
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^WGS  84 


^WGS  72  ^^*1* 


^GS  84  ”  ^GS  72 
^WGS  8  4  ”  ^WGS  72 


(7-1) 


As  indicated  from  Table  7.1,  when  proceeding  directly  from 
WGS  72  coordinates  to  obtain  WGS  84  values,  the  WGS  84  coordinates  will 
differ  from  the  WGS  72  coordinates  due  to  a  shift  in  the  coordinate 
system  origin,  a  change  in  the  longitude  reference,  a  scale  change 
(treated  through  Ar)  ,  and  changes  in  the  size  and  shape  of  the 
ellipsoid.  In  addition,  it  is  important  to  be  aware  that  A(j),  AX,  Ah 
calculated  using  Table  7.1  do  not  reflect  the  effect  of 
differences  between  the  WGS  72  and  WGS  84  EGMs  and  Geoids.  The 
following  cases  are  important  to  note: 

a.  Table  7.1  equations  are  to  be  used  for  direct  transformation 
of  Doppler-derived  WGS  72  coordinates.  These  transformed 
coordinates  should  agree  to  within  approximately  ±2  meters  with 
the  directly  surveyed  WGS  84  coordinates  using  TRANSIT  or  GPS 
point  positioning. 

b.  Table  7.1  should  not  be  used  for  satellite  local  geodetic 
stations  whose  WGS  72  coordinates  were  determined  using  datum 
shifts  from  [4],  The  preferred  approach  is  to  transform  such  WGS 
72  coordinates,  using  datum  shifts  from  [4],  back  to  their 
respective  local  datums,  and  then  transform  the  local  datum 
coordinates  to  WGS  84  using  Appendices  B  and  C. 

Table  7.1  should  be  used  only  when  no  other  approach  is 
applicable . 
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7 . 3  Local  Geodfitic . Pat.umttQrWGS -S.^-I^.tum-'IcaiiLS.fQX'mations, 


7.3.1  General 

Most  WGS  84  coordinates  needed  for  applications  and  DoD 
operations  in  Mapping,  Charting,  and  Geodesy  (MC&G)  will  be  obtained 
from  a  Local  Geodetic  Datum-to-WGS  84  Datum  Transformation.  This 
transformation  can  be  performed  either  in  curvilinear  (geodetic) 
coordinates : 


^WGS  84  “  ‘l^Local  ■*"  ^4* 

^GS  84  “  ^LocaJ  "*■  (7-2) 

^WGS  84  “  ^Local 

or,  in  rectangular  coordinates  [14]: 


X 

X 

Ax 

AS 

CO 

-¥ 

X  -  Xo 

Y 

Y 

+ 

Ay 

+ 

-CO 

AS 

F, 

o 

1 

(7-3) 

Z 

WGS  8  4 

Z 

Local 

Az 

¥ 

-e 

As 

Z  -  Zo 

Local 

where  AS  and  (e,  Yt  represent  changes  in  local  geodetic  datum  scale 
and  reference  frame  orientation,  respectively,  and  (Xq,  Yq,  Zq)  are  the 

coordinates  of  the  "initial"  (defining)  point  of  the  local  geodetic 
datum. 


There  are  several  datum  transformation  formulas  for 
accomplishing  the  preceding.  The  most  common  techniques  are,  in  the 
curvilinear  case,  the  Standard  MolodensJcy,  and  in  the  rectangular  case, 
the  3-,  4-,  or  7-parameter  transformations  depending  on  the 
availability  (and/or  reliability)  of  the  transformation  parameters. 
It  may  be  noted  that  the  3-parameter  rectangular  case  is  embedded 
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mathematically  in  the  Molodensky  Formulas  to  eliminate  the  conversion 
from  geodetic  to  rectangular  coordinates. 

In  addition,  the  curvilinear  and  rectangular  coordinate 
datum  transformations  can  be  accomplished  using  a  Multiple 
Regression  Equation  (MRE)  technique  which  accounts  for  the  non-linear 
distortion  in  the  local  geodetic  datum  [15] .  The  above  methods  are 
discussed  separately  in  [1] .  Only  the  Standard  Molodensky  Formula  and 
MRE  technique  are  discussed  here. 

7.3.2  The  Standard  Molodensky  Datum  Transformation  Formulas 


The  Standard  Molodensky  Datum  Transformation  Formulas 
[4]  [16],  along  with  definitions  of  the  terms,  are  listed  in  Table  7.2. 
As  the  Molodensky  Formulas  do  not  provide  satisfactory  results  near  the 
poles,  the  following  three-step  transformation  is  recommended: 


Local 


Kh) 


WGS  84 


(x,y. 


Z ) 

' Local 


I  (AX,AY,AZ) 


84 


Appendix  A  lists  the  reference  ellipsoid  names  and 
parameters  (semimajor  axis  and  flattening)  for  local  datums  currently 
tied  to  WGS  84  and  used  for  generating  datum  transformations. 

Appendix  B  contains  transformation  parameters  for  the 
geodetic  datums/systems  which  have  been  generated  from  satellite  ties 
to  the  respective  geodetic  control.  Due  to  the  errors  and  distortion 
that  affect  most  local  geodetic  datums,  use  of  mean  datum  shifts  (Ax, 
Ay,  AZ)  in  the  Standard  Molodensky  Datum  Transformation  Formulas  may 
produce  results  with  poor  quality  of  "fit".  Improved  fit  between  the 
local  datum  and  WGS  84  may  result  only  with  better  and  more  dense  ties 
with  local  or  regional  control  points. 
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Datum  transformation  shifts  derived  from  non-satellite 
information  are  available  in  Appendix  C. 

DMA-developed  local  geodetic  datum  geoid  heights  were 
used  in  forming  the  Local  Geodetic  Datum-to-WGS  84  Datum  Shifts  [17]  . 
An  example  of  such  a  geoid  for  NAD  27  is  included  in  contour  chart  form 
(Figure  7.1)  for  the  Contiguous  United  States  (CONUS). 

7.3.3  Datum  Transformation  Multiple  Regression  Equations 

The  development  of  Local  Geodetic  Datum-to-WGS  84  Datum 
Transformation  Multiple  Regression  Equations  [15]  was  initiated  to 
obtain  better  fits  over  continental  size  land  areas  than  could  be 
achieved  using  the  Standard  Molodensky  Formula  with  datura  shifts  (AX, 
Ay,  AZ) . 


For  A(|),  the  general  form  of  the  Multiple  Regression 
Equation  is  (also  see  [1] )  : 

A(i)  =  Aq  +  A^U  +  A2V  +  A3U2  +  A4UV  -I-  AsV=^  +  .  .  .  +  AggU^V®  (7-4) 


where 

Aq  =  constant 

A3,  Aj,  ...  =  coefficients  determined  in  the  development 

U  =  k  (  (|)  -  (|)^  )  =  normalized  geodetic  latitude  of  the 
computation  point 

V  =  k  (  A,  -  )  =  normalized  geodetic  longitude  of  the 

computation  point 

k  =  scale  factor,  and  degree-to-raa Lan  conversion 
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({),  X  =  local  geodetic  latitude  and  local  geodetic  longitude 
(in  degrees),  respectively,  of  the  computation  point 

<j)„,  X.jn  =  mid-latitude  and  mid-longitude  values,  respectively, 
of  the  local  geodetic  datum  area  (in  degrees) . 

Similiar  equations  are  obtained  for  and  Ah  by  replacing  A(j)  in  the 
left  portion  of  Equation  (7-4)  by  AX  and  Ah,  respectively. 

Local  Geodetic  Datum-to-WGS  84  Datum  Transformation 
Multiple  Regression  Equations  for  seven  major  continental  size  datums, 
covering  contiguous  continental  size  land  areas  with  large  distortion, 
are  provided  in  Appendix  D.  The  main  advantage  of  MRE '  s  lies  in 
modeling  of  distortion  for  better  fit  in  geodetic  applications. 
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Formulas  and  Parameters 
ransform  WGS  72  Coordinates 
to  WGS  84  Coordinates 
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Table  7 . 2 


Standard  Molodensky  Datum  Transformation  Formulas* 

-  Local  Geodetic  Datum  to  WGS  84  - 

1 .  The  Standand  Molodensky  Formulas 

A(j)"  =  {-  Ax  sin  (j)  cos  X  -  Ay  sin  (j)  sin  X,  +  AZ  cos  (j) 

+  Aa  (RNe^  sin  (j)  cos  (()) /a  +  Af  [R^, (a/b)  +  R^(b/a)]  sin  (|)  cos  (|)} 
•  [  (R-m  +  h)  sin  l"] 

AX"  =  [-  Ax  sin  X  +  Ay  cos  X]  •  [  (R^j  +  h)  cos  (j)  sin  l"]"^ 

Ah  =  Ax  cos  4>  cos  X  +  Ay  cos  (|)  sin  X  +  Az  sin  (|) 

-  Aa  (a/R,^)  +  Af  (b/a)  R|g  sin^(|) 


2 .  Definition  of  Terms  in  the  Molodensky  Formulas 


(j),  X,  h 


geodetic  coordinates  (old  ellipsoid) 


(j)  =  geodetic  latitude.  The  angle  between  the  plane  of  the 
geodetic  equator  and  the  ellipsoidal  normal  at  a  point 
(measured  positive  north  from  the  geodetic  equator, 
negative  south)  . 

X  =  geodetic  longitude.  The  angle  betvs’een  the  plane  of  the 
Zero  Meridian  and  the  plane  of  the  geodetic  meridian  of 
the  point  (measured  in  the  plane  of  the  geodetic 
equator,  positive  from  0°  to  180°  E,  and  negative 
from  0°  to  18i  ■’  W)  . 

h  =  geodetic  height  (ellipsoidal  height) .  The  distance  of  a 
point  from  the  ellipsoid  measured  from  the  surface  of 
the  ellipsoid  along  the  ellipsoidal  normal  to  the 
point . 

h  ~  N  +■  H 


N  =  ellipsoid  to  geoid  separation.  The  distance  of  the 
geoid  above  (+N)  or  below  (-N)  the  ellipsoid. 

H  =  distance  of  a  point  from  the  geoid  (or,  approximately, 
elevation  of  the  point  above/below  mean  sea  level); 
positive  above  mean  sea  level,  negative  below  mean  sea 
level . 


*  Not  to  be  used  between  89°  Latitude  and  the  pole  (see  Section  7.3.2) . 
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Table  7.2  (Cont'd) 


Standard  Molodensky  Datum  Transformation  Formulas 
-  Local  Geodetic  Datum  to  WGS  84  - 

A(j),  AX,  Ah  =  corrections  to  transform  local  geodetic  datum  coordi¬ 
nates  to  WGS  84  ([),  X,  h  values.  The  units  of  A(t)  and  AX. 
are  arc  seconds  (");  the  units  of  Ah  are  meters  (m) . 

NOTE:  AS  "h's"  ARE  NOT  AVAILABLE  FOR  LOCAL  GEODETIC 
DATUMS,  THE  Ah  CORRECTION  WILL  NOT  BE  APPLICABLE 
WHEN  TRANSFORMING  WGS84 . 

Ax,  Ay,  AZ  =  shifts  between  centers  of  the  local  geodetic  datum  and 
WGS  84  Ellipsoids;  corrections  to  transform  local 
geodetic  system-related  rectangular  coordinates  (X,  Y, 
Z)  to  WGS  84-related  X,  Y,  Z  values. 

a  “  semimajor  axis  of  the  local  geodetic  datum  ellipsoid. 

b  =  semiminor  axis  of  the  local  geodetic  datum  ellipsoid, 
b/a  =  1  -  f 

f  =  flattening  of  the  local  geodetic  datum  ellipsoid. 

Aa,  Af  =  differences  between  the  semimajor  axis  and  flattening 
of  the  local  geodetic  datum  ellipsoid  and  the  WGS  84 
Ellipsoid,  respectively  (WGS  84  minus  Local)  . 

e  =  first  eccentricity. 

e2  =  2f  -  f2 

=  radius  of  curvature  in  the  prime  vertical. 

=  a/ (l-e^sin^cj)) 

R,^  =  radius  of  curvature  in  the  meridian. 

Rm  =  a  (1-e^)  /  (l-e^sin^(})) 

NOTE :  All  A-quantities  are  formed  by  subtracting  local  geodetic  datum 

ellipsoid  values  from  WGS  84  Ellipsoid  values. 
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8.  ACCURACY  OF  WGS  84  COORDINATES 


The  accuracy  of  the  WGS  84  coordinates  of  a  site  is  significantly 
influenced  by  the  method  used  to  determine  the  coordinates.  Depending 
on  the  data  available,  the  WGS  84  coordinates  of  a  site  can  be 
determined: 

-  Directly  in  WGS  84  via  a  satellite  point  positioning  solution 
using  ground-based  Doppler  or  GPS  satellite  tracking  data  and 
broadcast  or  precise  satellite  ephemerides. 

-  By  a  WGS  72-to-WGS  84  Coordinate  Transformation. 

-  By  a  Local  Geodetic  Datum-to-WGS  84  Datum  Transformation. 

However,  the  situation  is  even  more  complicated  since  there  are 
several  techniques  for  accomplishing  a  Local  Geodetic  Datum-to-WGS  84 
Datum  Transformation.  In  addition,  the  accuracy  of  the  WGS  84 
coordinates  of  a  site  is  different  depending  on  whether  it  is  a 
satellite  station  or  a  non-satellite  geodetic  network  station,  or 
whether  the  WGS  84  coordinates  wej  e  determined  by  a  receiver  operated 
in  a  dynamic  or  static  mode. 

The  accuracy  (one  sigma)  of  WGS  84  coordinates  directly 
determined  in  WGS  84  by  Doppler  or  GPS  Satellite  Point  Positioning, 
their  respective  precise  ephemerides  and  ground-based  satellite 
tracking  data  acquired  in  the  static  mode,  is  in  geodetic  latitude, 
geodetic  longitude,  and  geodetic  height: 

“  Ox  =  ±1  m  (8-1) 

CTh  “  ±1  to  ±2  m.  (8-2) 
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The  Doppler  stations  used  in  the  development  of  WGS  84  were 
surveyed  prior  to  1  January  1987  in  the  NSWC  9Z-2  system.  As  such,  the 
indirectly  obtained  WGS  84  coordinates  (through  corrections  of  biases 
given  in  Chapter  2)  for  these  Doppler  stations  have  been  established  at 
lower  accuracies  compared  to  Doppler  stations  directly  surveyed  in  the 
WGS  84  reference  frame.  Thus,  the  absolute  accuracy  (one  sigma)  of 
these  Doppler  station  WGS  84  coordinates  was  assumed  to  be: 

“  ±2  m  (8-3) 

Oh  ”  ±2  to  ±3  m.  (8-4) 

The  WGS  84  coordinates  for  1591  Doppler  stations,  surveyed  up  to 
31  December  1986  and  used  as  an  integral  part  of  the  WGS  84 
development,  and  additional  stations,  surveyed  after  1  January  1987, 
have  been  used  to  develop  Local  Geodetic  Datum-to-WGS  84  Datum 
Transformations  (Chapter  7)  . 

The  WGS  84  coordinate  accuracies  in  the  two  paragraphs  immediately 
above  are  absolute  accuracies  in  that  they  incorporate  not  only  the 
"observational"  or  solution  error,  but  the  errors  associated  with 
placing  the  origin  of  the  WGS  84  Coordinate  System  at  the  earth's 
center  of  mass  and  determining  the  correct  scale  for  the  WGS  84 
Coordinate  System.  The  error  estimates  do  not  include  the  uncertainty 
associated  with  the  attempt  to  bring  the  WGS  84  zero  meridian  into 
coincidence  with  the  BIH-defined  Zero  Meridian  for  the  epoch  1984.0. 
This  is  not  necessary  since  the  location  of  the  WGS  84  longitude 
reference  or  zero  meridian  is  arbitrary.  These  absolute  accuracy  values 
should  not  be  confused  v;ith  the  sub-meter  precision: 

-  Of  a  Doppler  or  GPS  coordinate  solution  (the  "observi  tional" 
error) . 

-  Of  a  Doppler  or  GPS  coordinate  solution  which  has  been  repeated 
independently  at  the  same  site. 
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The  WGS  84  coordinates  of  a  non-satellite  derived  local  geodetic 
network  station  will  be  less  accurate  than  the  WGS  84  coordinates  of  a 
Doppler  or  GPS  station.  This  is  due  to  the  distortions  and  surveying 
errors  present  in  local  geodetic  datum  networks,  the  lack  (in  general) 
of  a  sufficient  number  of  properly  placed  Doppler  or  GPS  stations 
colocated  with  local  geodetic  datum  stations  for  use  in  forming  the 
Local  Geodetic  Datum-to-WGS  84  Datum  Shifts,  and  the  uncertainty 
introduced  by  the  datum  transformation. 
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9.  CONCLUSIONS /SUMMARY 


World  Geodetic  System  1984  is  based  on  the  use  of  data, 
techniques,  and  technology  available  in  early  1984.  As  a  result,  WGS  84 
is  more  accurate  than  WGS  72  and  replaces  the  latter  as  the  geocentric 
system  officially  authorized  for  DoD  use. 

The  origin  and  orientation  of  the  WGS  84  Reference  Frame  are  more 
accurately  defined  than  they  were  for  WGS  72.  In  addition,  Doppler  and 
GPS-derived  Local  Geodetic  Datum-to-WGS  84  Datum  Shifts  are  more 
accurate  than  analogous  WGS  72  values,  and  are  available  for  many  more 
datums  for  WGS  84  as  compared  to  WGS  72.  Further,  the  WGS  84  EGM  and 
geoid  are  considerably  more  accurate  than  their  WGS  72  counterparts, 
and  minor  scale  errors  inherent  in  WGS  72  are  reduced  in  WGS  84.  These 
improvements  translate  into: 

-  More  accurate  maps  and  charts  of  scale  1:50,000  and  larger. 

-  More  accurate  geodetic  coordinates,  geoid  heights,  heights 
above  the  geoid  (approximately  mean  sea  level),  and  distances. 

-  An  improved  capability  for  satellite  orbit  determination  and 
prediction . 

-  The  capability  to  place  many  more  local  geodetic  datums  on  a 
world  geodetic  system,  and  do  it  more  accurately. 

The  latter  is  particularly  important  for  those  local  geodetic 
datums  affected  by  large  distortions.  Placement  of  such  local  datums  on 
WGS  84,  using  the  variable  datum  shifts  made  possible  by  a  well 
dispersed  set  of  Doppler  or  GPS  sites,  effectively  removes  these 
distortions . 
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The  value  of  WGS  84  will  become  increasingly  evident  in  the  early 
1990s  when  Navstar  GPS  will  be  fully  operational.  Since  the  reference 
system  for  Navstar  GPS  is  WGS  84,  high  quality  geocentric  coordinates 
can  be  provided  automatically  by  Navstar  GPS  User  Equipment.  For  those 
using  Navstar  GPS  but  still  utilizing  local  geodetic  datums  and 
products,  the  availability  of  the  more  accurate  WGS  84-to-Local 
Geodetic  Datum  Shifts  will  lead  to  an  improved  recovery  of  local 
coordinates.  Again,  the  value  of  having  all  MC&G  products  and 
navigational  activities  referenced  to  WGS  84  is  noted.  But  if  local 
geodetic  datums  are  in  use,  requiring  a  WGS  84-to-Local  Geodetic  Datum 
Transformation,  then  the  value  of  havj.ng  improved  datum  shifts  (made 
possible  by  a  well  dispersed  set  of  Doppler  or  GPS  sites  throughout  the 
region)  is  apparent. 

Efforts  have  been  initiated  to  enhance/refine  WGS  84  to  satisfy 
anticipated  future  requirements  for  MC&G  products  and  data  of  increased 
accuracy . 
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APPENDIX  A 


LIST  OF  REFERENCE  ELLIPSOID  NAMES  AND  PARAMETERS 
(USED  FOR  GENERATING  DATUM  TRANSFORMATIONS) 
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REFERENCE  ELLIPSOIDS 
FOR 

LOCAL  GEODETIC  DATUMS 


1 ,  GENERAL 

This  appendix  lists  the  reference  ellipsoids  and  their 
constants  (a,f)  associated  with  the  local  geodetic  datums  which  are 
tied  to  WGS  84  through  datum  transformation  constants  and/or  MRE's 
(Appendices  B,  C,  and  D) . 

2.  CONSTANT  CHARACTERSTICS 

In  Appendix  A.l,  the  list  of  ellipsoids  includes  a  new 
feature.  Some  of  the  reference  ellipsoids  have  more  than  one  semi¬ 
major  axis  (a)  associated  with  them.  These  different  values  of  axis 
(a)  vary  from  one  region  or  country  to  another  or  from  one  year  to 
another  within  the  same  region  or  country. 

A  typical  example  of  such  an  ellipsoid  is  Everest  whose 
semi-major  axis  (a)  was  originally  defined  in  yards.  Here,  changes 
in  the  yard  to  meter  conversion  ratio  over  the  years  have  resulted 
in  five  different  values  for  the  constant  (a),  as  identified  in 
Appendix  A.l. 

To  facilitate  correct  referencing,  a  standardized  two 
letter  code  is  also  included  to  identify  the  different  ellipsoids 
and/or  their  "versions"  pertaining  to  the  different  values  of  the 
semi-major  axis  (a). 
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Reference  Ellipsoid  Names  and  Constants 
Used  for  Datum  Transformations* 
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As  accepted  by  DMA. 

Thrcugh  adoption  of  a  new  yard  to  meter  conversion  factor  in  the  referenced  country. 
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APPENDIX  B 

DATUM  TRANSFORMATIONS  DERIVED 
USING  SATELLITE  TIES  TO  GEODETIC  DATUMS/SYSTEMS 
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DATUM  TRANSFORMATION  CONSTANTS 
-  GEODETIC  DATUMS /SYSTEMS  TO  WGS  84 
(THROUGH  SATELLITE  TIES) 


1 .  GENERAL 

This  appendix  provides  the  details  about  the  reference 
ellipsoids  (Appendix  A)  which  are  used  as  defining  parameters  for 
the  geodetic  datums  and  systems . 

There  are  99  local  geodetic  datums  which  are  currently 
related  to  WGS  84  through  satellite  ties. 

2.  LOCAL  DATUM  ELLIPSOIDS 

Appendix  B.l  lists,  alphabetically,  the  local  geodetic 
datums  and  the  Soviet  Geodetic  system  1985  (SGS  85)  with  their 
associated  ellipsoids.  Two  letter  ellipsoidal  codes  (Appendix  A) 
have  also  been  included  to  clearly  indicate  which  "version"  of  the 
ellipsoid  was  used  in  determining  the  transformation  constants. 

3 .  TRANSFORMATION  CONSTANTS 

Appendices  B.2  through  B.7  list  the  constants  for  local 
datums  for  continental  areas.  The  continents  and  the  local  geodetic 
datums  are  arranged  alphabetically. 

Appendices  E.8  through  B.IO  list  the  constants  for  local 
datums  which  fall  within  the  ocean  areas.  The  ocean  areas  and  the 
geodetic  datums  are  also  arranged  alphabetically. 


4 .  ERROR  ESTIMATES 

The  la  error  estimates  for  the  datum  transformation 
constants  (AX,Ay,AZ),  obtained  from  the  computed  solutions,  are 
also  tabulated.  These  estimates  do  not  include  the  errors  of  the 
common  control  station  coordinates  which  were  used  to  compute  the 
shift  constants. 

For  datums  having  four  or  less  common  control  stations, 
the  la  errors  for  shift  constants  are  non-computed  estimates. 

The  current  set  of  error  estimates  have  been  revaluated 
and  revised  after  careful  consideration  of  the  datum  transformation 
solutions  and  the  related  geodetic  information;  the  intent  has  been 
to  assign  estimates  as  realistic  as  possible. 
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Appendix  B.l 


Geodetic  Datums/Reference  Systems 
Related  to  World  Geodetic  System  1984 
(Through  Satellite  Ties) 


Local  Geodetic  Datum 


Adindan 

Afgooye 

Aiii  el  Abd  1970 
Anna  1  Astro  1965 
Antigua  Island  Astro  1943 

Arc  1950 
Arc  1960 

Ascension  Island  1958 
Astro  Beacon  "E"  1945 
Astro  DOS  71/4 

Astro  Tern  Island  (FRIG)  1961 
Astronomical  Station  1952 
Australian  Geodetic  1966 
Australian  Geodetic  1984 
Ayabelle  Lighthouse 

Bellevue  (IGN) 

Bermuda  1957 
Bissau 

Bogota  Observatory 
Campo  Inchauspe 

Canton  Astro  1966 
Cape 

Cape  Canaveral 
Carthage 

Chatham  Island  Astro  1971 

Chua  Astro 
Corrego  Alegre 
Dabola 

Djakarta  (Batavia) 

DOS  1968 


Associated* 


Reference  Ellipsoid 

Code 

Clarke  1880 

CD 

Krassovsky  1940 

KA 

International 

1924 

IN 

Australian  National 

AN 

Clarke  1880 

CD 

Clarke  1880 

CD 

Clarke  1880 

CD 

International 

1924 

IN 

International 

1924 

IN 

International 

1924 

IN 

International 

1924 

IN 

International 

1924 

IN 

Australian  National 

AN 

Australian  National 

AN 

Clarke  1880 

CD 

International 

1924 

IN 

Clarke  1866 

CC 

International 

1924 

IN 

International 

1924 

IN 

International 

1924 

IN 

International 

1924 

IN 

Clarke  1880 

CD 

Clarke  1866 

CC 

Clarke  1880 

CD 

International 

1924 

IN 

Internationax 

1924 

IN 

International 

1924 

IN 

Clarke  1880 

CD 

Bessel  1841 

BR 

International 

1924 

IN 

*  See  Appendix  A.l  for  associated  constants  a,f. 
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Appendix  B.l  (Cont'd) 


Geodetic  Datums/Reference  Systems 
Related  to  World  Geodetic  System  1984 
(Through  Satellite  Ties) 


Local  Geodetic  Datum 

Associated* 

Reference  Ellipsoid 

Code 

Easter  Island  1967 

International  1924 

IN 

European  1950 

International  1924 

IN 

European  1979 

International  1924 

IN 

Fort  Thomas  1955 

Clarke  1880 

CD 

Gan  1970 

International  1924 

IN 

Geodetic  Datum  1949 

International  1924 

IN 

Graciosa  Base  SW  1948 

International  1924 

IN 

Guam  1963 

Clarke  1866 

CC 

GUX  1  Astro 

International  1924 

IN 

Hjorsey  1955 

International  1924 

IN 

Hong  Kong  1963 

International  1924 

IN 

Hu-Tzu-Shan 

International  1924 

IN 

Indian 

Everest 

EA/EC* 

Indian  1954 

Everest 

EA 

Indian  1975 

Everest 

EA 

Ireland  1965 

Modified  Airy 

AM 

ISTS  061  Astro  1968 

International  1924 

IN 

ISTS  073  Astro  1969 

International  1924 

IN 

Johnston  Island  1961 

International  1924 

IN 

Kandawala 

Everest 

EA 

Kerguelen  Island  1949 

International  1924 

IN 

Kertau  1948 

Everest 

EE 

Kusaie  Astro  1951 

International  1924 

IN 

L.  C.  5  Astro  1961 

Clarke  1866 

CC 

Leigon 

Clarke  1880 

CD 

Liberia  1964 

Clarke  1880 

CD 

Luzon 

Clarke  1866 

CC 

Mahe  1971 

Clarke  1880 

CD 

Massawa 

Bessel  1841 

BR 

Merchich 

Clarke  1880 

CD 

*  See  Appendix  A.l  for  associated  constants  a,f. 

**  Due  to  different  semi-major  axes.  See  Appendix  A.l. 
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Appendix  B.l  (Cont'd) 


Geodetic  Datums/Reference  Systems 
Related  to  World  Geodetic  System  1984 
(Through  Satellite  Ties) 


Local  Geodetic  Datum  Associated* 


Reference  Ellipsoid 

Midway  Astro  1961 

International 

1924 

IN 

Minna 

Clarke  1880 

CD 

Montserrat  Island  Astro  1958 

Clarke  1880 

CO 

M'Poraloko 

Clarke  1880 

CD 

Nahrwan 

Clarke  1880 

CD 

Naparima,  BWI 

International 

1924 

IN 

North  American  1927 

Clarke  1866 

CC 

North  American  1983 

GRS  80*** 

RF 

Observatorio  Meteorologico  1939 

International 

1924 

IN 

Old  Egyptian  1907 

Helmert  1906 

HE 

Old  Hawaiian 

Clarke  1866 

CC 

Oman 

Clarke  1880 

CD 

Ordnance  Survey  of  Great  Britain  1936 

Airy  1830 

AA 

Pico  de  las  Nieves 

International 

1924 

IN 

Pitcairn  Astro  1967 

International 

1924 

IN 

Point  58 

Clarke  1880 

CD 

Pointe  Noire  1948 

Clarke  1880 

CD 

Porto  Santo  1936 

International 

1924 

IN 

Provisional  South  American  1956 

International 

1924 

IN 

Provisional  South  Chilean  1963**** 

International 

1924 

IN 

Puerto  Rico 

Clarke  1866 

CC 

Qatar  National 

International 

1924 

IN 

Qornoq 

International 

1924 

IN 

Reunion 

International 

1924 

IN 

Rome  1940 

International 

1924 

IN 

Santo  (DOS)  1965 

International 

1924 

IN 

Sao  Braz 

International 

1924 

IN 

Sapper  Hill  1943 

International 

1924 

IN 

Schwarzeck 

Bessel  1841 

BN 

Selvagem  Grande  1938 

International 

1924 

IN 

*  See  Appendix  A.l  for  associated  constants  a,f. 
***  Geodetic  Reference  System  1980 
****Also  known  as  Hito  XVIII  1963 
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Appendix  B.l  (Cont'd) 


Geodetic  Datums/Reference  Systems 
Related  to  World  Geodetic  System  1984 
(Through  Satellite  Ties) 


Local  Geodetic  Datum 


South  American  1969 
South  Asia 
Timbalai  1948 
Tokyo 

Tristan  Astro  1968 

Viti  Levu  1916 
Wake-Eniwetok  1960 
Wake  Island  Astro  1952 
Zander! j 


See  Appendix  A.l  ^or  associated 


Associated* 

Reference  Ellipsoid  Code 

South  American  1969  SA 
Modified  Fischer  1960  FA 
Everest  EB 
Bessel  1841  BR 
International  1924  IN 

Clarke  1880  CD 
Hough  1960  HO 
International  1924  IN 
International  1924  IN 
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Datum  1950  coordinates  developed  from  Ordnance  Survey  of  Great  Britain  (OSGB)  Scientific 
1980  (SN  80)  coordinates. 
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Appendix  C.l 


Local  Geodetic  Datums 
Related  to  World  Geodetic  System  1984 
(Through  non-Satellite  Ties) 


Local  Geodetic  Datum 

Assoc^iated* 
Reference  Ellipso.id 

Code 

Bukit  Rimpah 

Bessel  1841 

BR 

Camp  Area  Astro 

International  1924 

IN 

European  1950 

International  1924 

IN 

Gunung  Segara 

Bessel  1841 

BR 

Herat  North 

International  1924 

IN 

Tananarive  Observatory  1925 

International  1924 

IN 

Yacare 

International  1924 

IN 

*  See  Appendix  A.l  for  associated  constants  a,f. 
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DATUM  TRANSFORMATION  CONSTANTS 
LOCAL  GEODETIC  DATUMS  TO  WGS  84 
(THROUGH  NON-SATELLITE  TIES) 


1 .  GENERAL 

This  appendix  provides  the  details  about  the  reference 
ellipsoids  (Appendix  A)  used  as  defining  parameters  for  the  local 
geodetic  datums  which  are  related  to  WGS  84  through  non-satellite 
ties  to  the  local  control. 

There  are  six  such  local  geodetic  datums,  and  one  special 
area  under  the  European  Datum  1950  (ED  50). 

2.  LOCAL  DATUM  ELLIPSOIDS 

Appendix  C.l  lists  alphabetically  the  local  geodetic 
datums  and  their  associated  ellipsoids.  Two  letter  ellipsoidal 
codes  (Appendix  A)  have  also  been  included  to  clearly  indicate 
which  "version"  of  the  ellipsoid  has  been  used  to  determine  the 
transformation  constants. 

3 .  TRANSFORMATION  CONSTANTS 

Appendix  C.2  alphabetically  lists  the  local  geodetic 
datums  and  the  special  area  under  ED  50  with  the  associated  shift 
constants . 


4 .  ERROR  ESTIMATES 

The  error  estimates  are  not  available  for  the  datum 
transformation  constants  listed  in  the  Appendix  C.2. 
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FOR 
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MULTIPLE  REGRESSION  EQUATIONS 


1 .  GENERAL 

This  appendix  provides  the  Multiple  Regression  Equations 
(MRE's)  parameters  for  continental  size  datums  and  for  contiguous 
large  land  areas  (Table  D.l). 

Table  D.l 

DATUMS  WITH  MULTIPLE  REGRESSION  EQUATIONS 


DATUM  NAME 

AREA  COVERED 

Australian  Geodetic 

1966 

Australian  Mainland 

Australian  Geodetic 

1984 

Australian  Mainland 

Campo  Inchauspe 

Argentina 

Corrego  Alegre 

Brazil 

European  1950 

Western  Europe 

Austria,  Denmark,  France, 
W.  Germany*,  The  Nether¬ 
lands,  and  Switzerland. 

North  American  1927 

CONUS 

Canadian  Mainland 

South  American  1969 

South  American  Mainland 
Argentina,  Bolivia, 

Brazil,  Chile,  Colombia, 
Ecuador,  Guyana,  Peru, 
Paraguay,  Uruguay,  and 
Venezuela . 

*  Prior  to  October  1990. 

2 .  APPLICATIONS 

The  coverage  area  for  MRE's  application  are  defined  in 
detail  in  Appendices  D.l  through  D.6.  MRE's  coverage  area  should 
never  be  extrapolated  and  are  not  to  be  used  over  islands  and/or 
isolated  land  areas. 

The  main  advantage  of  MRE’s  lies  in  their  modeling  of 
distortions  for  datums,  which  cover  continental  size  land  areas,  to 
obtain  better  transformation  fit  in  geodetic  applications. 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

Australian  Geodetic  Datum  1966  (AUA)  to  WGS  84 


Area  of  Applicability  :  Australian  Mainland  (excluding  Tasmania) 


MRE  coefficients  for  4>  and  A,  are  : 

A(|)"  =  5.19238  +  0.12666  U  +  0.52309  V  -  0.42069 

-  0.39326  UV  +  0.93484  uV  +  0.44249  UV^  -  0.30074  UV^ 

+  1.00092  -  0.07565  V®  -  1.42988  U®  -  16.06639  uV® 

+  0.07428  V®  +  0.24256  UV®  +  38.27946  uV 

-  62.06403  uV  +  89.19184  U®V® 


AA,"  =  4.69250  -  0.87138  U  -  0.50104  V  +  0.12678  UV 

-  0.23076  V"*  -  0.61098  uV  -  0.38064  V^  +  2.89189  U® 

+  5.26013  uV  -  2.97897  U®  +  5.43221  uV 

-  3.40748  uV  +  0.07772  V®  +  1.08514  U®V  +  0.71516  UV® 

+  0.20185  V®  +  5.18012  uV®  -  1.72907  uV 

-  1.24329  uV 


where  :  U  =  K  (4>  +  27°);  V  =  K  (A,  -  134°);  K  =  0.05235988 


NOTE  ; 


Input  <j)  as  (-)  from  90°S  to  0°N  in  degrees. 
Input  A  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  ±  2.0  m 


Test  Case  : 


AUA 


Shift 


WGS  84 


<{)=  (-)  17°00' 32.78"S  A<j)=  5.48" 

A=  144°11' 37.25"E  AA=  3.92" 


(J)=  (-)  17°00'27.30"S 
A=  144°11'41. 17"E 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

Australian  Geodetic  Datum  1984  (AUG)  to  WGS  84 


Area  of  Applicability  :  Australian  Mainland  (excluding  Tasmania) 


MRE  coefficients  for  (|)  and  X  are  : 


A(l)"  =  5.20604  +  0.25225  U  +  0.58528  V  -  0.41584 

-  0.38620  UV  -  0.06820  +  0.38699  uV  +  0.07934  UV^ 
+  0.37714  -  0.52913  uV  +  0.38095  v’  +  0.68776  uV 

-  0.03785  V®  -  0.17891  U®  -  4.84581  uV  -  0.35777  v’ 

+  4.23859  uV 


AX"  =  4.67877  -  0.73036  U  -  0.57942  V  +  0.28840 

+  0.10194  -  0.27814  UV^  -  0.13598  +  0.34670  UV^ 

-  0.46107  V*  +  1.29432  uV  +  0.17996  UV^  -  1.13008  U 

-  0.46832  U®  +  0.30676  V®  +  0.31948  U®  +  0.16735  V® 

-  1.19443  uV 


Where  :  U  =  K  (<{)  +  27°);  V  =  K  (X  -  134°);  K  =  0.05235988 


NOTE  :  Input  <J)  as  (-)  from  90°S  to  0°N  in  degrees. 

Input  X  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  i  2.0  m 


Test  Case  : 


AUG 


Shift 


WGS  84 


<J)=  (-)20°38'00.67"S  A(t>=  5.50" 
X=  144°24 '29.29"E  AX=  4.11" 


(J)=  (-)20°37’55.17"S 
X=  144°24'33.40"E 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

Campo  Inchauspe  Datum  (CAI)  to  WGS  84 


Area  of  Applicability  :  Argentina  (Continental  land  areas  only) 


MRE  coefficients  for  <|)  and  A,  are  : 

A4)"  =  1.67470  +  0.52924  U  -  0.17100  V  +  0.18962 

+  0.04216  UV  +  0.19709  UV^  -  0.22037  U*  -  0.15483  uV 

-  0.24506  UV*  -  0.05675  V®  +  0.06674  U®  +  0.01701  UV® 

-  0.00202  U’  +  0.08625  v’  -  0.00628  U®  +  0.00172  U®V* 

+  0.00036  uV 

AA,"  =  -  2.93117  +  0.18225  U  +  0.69396  V  -  0.04403  U^ 

+  0.07955  V^  +  1.48605  V®  -  0.00499  U*  -  0.02180  uV 

-  0.29575  uV  +  0.20377  UV*  -  2.47151  V®  +  0.09073  uV 
+  1.33556  v’  +  0.01575  uV  -  0.26842  V® 


Where  :  U  =  K  (  4»  +  35“);  V  =  K  (  A,  +  64°);  K  =  0.15707963 


NOTE  :  Input  <t>  as  (-)  from  90°S  to  0°N  in  degrees. 

Input  A  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  i  2.0  m 


Test  Case  ; 

CAI  Shift  WGS  84 

4>=  (-)29°47'45.68"S  A<J»=  1.95"  <J)=  (  -  )  29°47 ' 43 . 73"S 

A=  (-)58°07 ' 38.20"W  AA=-1.96"  A=  ( - ) 58°07 ' 40 . 16"W 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

Corrego  Alegre  Datum  (COA)  to  WGS  84 


Area  of  Applicability  :  Brazil  (Continental  land  areas  only) 


MRE  coefficients  for  <|)  and  X  are  : 

A<t)"  =  -  0.84315  +  0.74089  U  -  0.21968  V  -  0.98875  u" 

+  0.89883  UV  +  0.42853  +  2.73442  -  0.34750  uV 

+  4.69235  -  1.87277  +  11.06672  U^V 

-  46.24841  uV  -  0.92268  u’  -  14.26289  uV 

+  334.33740  -  15.68277  uV  -  2428.8586  uV 

Al"  =  -  1.46053  +  0.63715  U  +  2.24996  V  -  5.66052  UV 

+  2.22589  -  0.34504  U^  -  8.54151  uV  +  0.87138  U* 

+  43.40004  uV  +  4.35977  UV^  +  8.17101  uV 
+  16.24298  uV  +  19.96900  UV'‘  -  8.75655 

-  125.35753  U®V  -  127.41019  uV  -  0.61047  U® 

+  138.76072  U'v  +  122.04261  uV  -  51.86666  uV 
+  45.67574 


Where  :  U  =  K  (  4)  +  15°);  V  =  K  (  A,  +  50°);  K  =  0.05235988 

NOTE  :  Input  <j)  as  (-)  from  90°S  to  0°N  in  degrees. 

Input  A,  as  (-)  from  180°W  to  0°K  in  degrees. 


Quality  of  fit  =  i  2.0  m 


Test  Case  : 

COA  Shift  WGS  84 

<l>=  (-)20°29 ' 01 .02"S  A<1)=  -1.03"  4)=  (  -  )  20°29 ' 02 . 05"S 

X=  (-)54°47 ' 13 . 17"W  AA-  -2.10"  A=  ( - ) 54°47 ' 15 . 27"W 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

European  Datum  1950  (EUR)  to  WGS  84 


Area  of  Applicability  :  Western  Europe*  (Continental  contiguous 

land  areas  only) 


MRE  coefficients  for  <J)  and  A  are  : 

A<t)"  =  -  2.65261  +  2.06392  U  +  0.77921  V  +  0.26743 

+  0.10706  UV  +  0.76407  -  0.95430  uV  +  0.17197  U* 

+  1.04974  uV  -  0.22899  uV  -  0.05401  V®  -  0.78909  U® 

-  0.10572  uV  +  0.05283  UV®  +  0.02445  uV 

AA"  =  -  4.13447  '  1.50572  U  +  1.94075  V  -  1.37600  U^ 

+  1.98425  UV  +  0.30068  V^  -  2.31939  U®  -  1.70401  U* 

-  5.48711  UV®  4-  7.41956  U®  ~  1.61351  uV 

+  5.92923  UV'*  -  1.97974  V^  +  1.57701  U®  -  6.52522  U®V® 
+  16.85976  uV  -  1.79701  UV®  -  3.08344  U^ 

-  14.32516  U®V  +  4.49096  uV^  +  9.98750  U®V 
+  7.80215  uV  -  2.26917  uV'  +  0.16438  V® 

-  17.45428  uV  -  8.25844  U®V®  +  5.28734  U®V® 

+  8.87141  U®V’  -  3.48015  U®V^  +  0.71041  uV 


Where  :  U  =  K  (  <|)  -  52°);  V  =  K  (  A  -  10°);  K  =  0.05235988 

NOTE  ;  Input  (|)  as  (-)  from  90°S  to  0°N  in  degrees. 

Input  A  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  +  2.0  m 


Test  Case  : 

EUR  Shift  WGS  84 

4»=  46°41 ' 42 . 89"N  A(j)=  -3.08"  <|>=  46 °4 1 ' 39 . 81  "N 

A=  13°54 ' 54 . 09"E  AA=  -3.49"  A=  13°54 ' 50 . 60"E 


*  See  Table  D.l  (Page  D-3)  for  the  list  of  countries 
covered  by  the  above  set  of  MRE's. 


Multiple  Regression  Equations  (MRE's) 
for  Transforming 

North  American  Datum  1927  (NAS)  to  WGS  84 

Area  of  Applicability  :  Canada  (Continental  contiguous  land 

areas  only) 


MRE  coefficients  for  4)  and  X  are  : 


A<|)" 


AX" 


0.79395  +  2.29199  U  +  0.27589  V  -  1.76644  u" 

+  0.47743  UV  +  0.08421  -  6.03894  -  3.55747 

-  1.81118  UV"*  -  0.20307  v'  +  7.75815  -  3.1017  uV 

+  3.58363  uV  -  1.31086  UV^  -  0.45916  +  14.27239  U® 

+  3.28815  \j\  +  1.35742  uV  +  1.75323  UV^  +  0.44999  V® 

-  19.02041  uV^  -  1.01631  uV  +  1.47331  UV® 

+  0.15181  V®  +  0.41614  uV  -  0.80920  UV®  -  0.18177  V^ 

+  5.19854  uV  -  0.48837  UV"  -  0.01473  V®  -  2.26448  U® 

-  0.46457  uV  +  0.11259  UV®  +  0.02067  V® 

+  47.64961  U®V^  +  0.04828  UV®  +  36.38963  U®V^ 

+  0.06991  uV  +  0.08456  uV  +  0.09113  uV 

+  5.93797  uV  -  2.36261  uV  +  0.09575  U®V® 

-  1.36099  +  3.61796  V  -  3.97703  U^  +  3.09705  UV 

-  1.15866  V^  -  13.28954  U®  -  3.15795  U®V  +  0.68405  UV® 

-  0.50303  V®  -  8.81200  U®V  -  2.17587  U®V®  -  1.49513  UV 

+  0.84700  V*  +  31.42448  U®  -  14.67474  U®V® 

+  0.65640  UV*  +  17.55842  U®  +  6.87058  uV  -  0.21565  V® 

+  62.18139  U®V®  +  1.78687  U®V*  +  2.74517  U®V® 

-  0.30085  UV®  +  0.04600  v’  +  63.52702  U®V® 

+  7.83682  U®V®  +  9.5944*  U®V®  +  0.01480  V®  + 

+  10.51228  uV®  -  1,42398  U®v"  -  0.00834  V® 

+  5.23485  uV  -  3.18129  U®v’  +  8.45704  U®V® 

-  2.29333  uV  +  0.14465  U®V®  +  0.29701  U®V® 

+  0.17655  uV 


Where  :  U  =  K  (  <J)  -  60°);  V  =  K  (  X  +  100°);  K  =  0.05235988 


NOTE  ; 


Input  <t)  as  (-)  from  90°S  to  0°N  in  degrees. 
Input  X  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  +  2,0  m 


Test  Case  ; 

NAS 


Shift 


WGS  84 


4)=  54°26'08.67"N 

X=  (-)110°17’02.41"W 


A<|)=  0.29"  4)=  54°26'08.96"N 

AX=  -3.16"  X=  (-)110°17'05.57"W 
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Multiple  Regression  Equations  (MRE's) 
for  Transforming 

North  American  Datum  1927  (NAS)  to  WGS  84 


Area  of  Applicability  :  USA  (Continental  contiguous  land  areas 

only;  excluding  Alaska  and  Islands) 


MRE  coefficients  for  4)  and  1  are  ; 


A4>" 


AA" 


0.16984  - 

-  4.57801 
+  0.12415 
+  0.73357 
+  0.07439 

-  1.30575 


0.76173  U  +  0.09585  V  +  1.09919  u" 

-  1.13239  uV  +  0.49831  -  0.98399  uV 

UV^  +  0.11450  +  27.05396  U®  +  2.03449  uV 
uV-  0.37548  -  0.14197  V®  -  59.96555  u" 

-  4.76082  U®  +  0.03385  V®  +  49.04320  U® 
U®V®  -  0.07653  uV+  0.08646  uV 


-  0.88437  +  2.05061  V  +  0.26361  U  -  0.76804  UV 

I  /\  ^2  .|  Onn*1tA  TT^  A  rrAI/TA  TT^tT  1  ACTAC' 


-  0.49211 

uV' 

+  2.17204 

UV®  - 

0.06004 

V^  +  0.30139 

uV 

+  1.88585 

uv^ 

-  0.81162 

UV®  - 

0.05183 

V®  -  0.96723 

UV® 

-  0.12948 

uV® 

+  3.41827 

U®  - 

0.44507 

U°V  +  0.18882 

UV 

0.01444  V®  +  0.04794  UV®  -  0.59013  U®V® 


Where : 


U  =  K  (  <|)  -  37°);  V  =  K  (  A  +  95°);  K  =  0.05235988 


NOTE  : 


Input  4)  as  (-)  from  90°S  to  0°N  in  degrees. 
Input  A  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  ±  2.0  m 


Test  Case  : 


NAS 


Shift 


WGS  84 


4)=  34°47 '08.83"N  A4>=  0.36" 

A=  (-)  86°34'52.18"W  AA=  0.08" 


4>=  34°47'09.19"N 

A=  (-)  86°34'52.10"W 


D-11 


Multiple  Regression  Equations  (MRE's) 
for  Transforming 

South  American  Datum  1969  (SAN)  to  WGS  84 


Area  of  Applicability  :  South  America  (Continental  contiguous 

]and  areas  only) 


MRE  coefficients  for  4)  and  X  are  : 

A4)”  =  -  1.67504  -  0.05209  U  +  0.25158  V  +  1.10149 

+  0.24913  UV  -  1.00937  uV  -  0.74977  -  1.54090  U* 

+  0.14474  +  0.47866  +  0.36278  uV  -  1.29942  UV^ 

+  0.30410  +  0.87669  U®  -  0.27950  uV  -  0.46367  u’ 

+  4.31466  UV^  +  2.09523  uV  +  0.85556  UV®  -  0.17897  U® 

-  0.57205  Uv"  +  0.12327  U®  -  0.85033  U®V®  -  4.86117  uV 
+  0.06085  uV  -  0.21518  uV®  +  0.31053  uV 

-  0. 09228  U®V®  -  0.22996  uV  +  0.58774  U®V® 

+  0.87562  uV  +  0.39001  U®V®  -  0.81697  uV 

AA,"  =  -  1.77967  +  0.40405  U  +  0.50268  V  -  0.05387  u" 

-  0.12837  UV  -  0.54687  uV  -  0.17056  V®  -  0.14400  uV 

+  0.11351  U®V  -  0.62692  uV®  -  0.01750  U®  +  1.18616  uV 
+  0.01305  U®  +  1.01360  uV  -  0.29059  U®V® 

+  5.12370  U®V®  -  5.09561  uV  -  5.27168  U®v' 

+  4.04265  uV  -  1.62710  U®v’  +  1.68899  uV 

+  2.07213  U®V®  -  1.76074  uV 


Where:  U  =  K  (  4)  +  20°);  V  =  K  (  A.  +  60°);  K  =  0.05235988 


NOTE  :  Input  4)  as  (-)  from  90°S  to  0°N  in  degrees. 

Input  A,  as  (-)  from  180°W  to  0°E  in  degrees. 


Quality  of  fit  =  +  2.0  m 


Test  Cc\se  : 

SAN  Shift  WGS  84 

4)-  (-)  31°56'33.95"S  A4>=  -1.36"  4)=  (-)  31°56 ' 35 . 31"S 

A.=  (-)  65°06 ' 18.66"W  AA=  -2.16"  A,=  (-)  65°06 ' 20 . 82  "W 
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APPENDIX  E 
-ACRONYMS- 


AGD  66 

— 

Australian  Geodetic  Datum  1966 

AGD  84 

- 

Australian  Geodetic  Datum  1984 

BIH 

= 

Bureau  International  de  I'Heure 

BTS 

ts 

BIH  Terrestrial  System 

CEP 

= 

Celestial  Ephemeris  Pole 

CIS 

Conventional  Inertial  System 

CONUS 

SB 

Contiguous  United  States 

CTP 

= 

Conventional  Terrestrial  Pole 

CTS 

SS 

Conventional  Terrestrial  System 

DMA 

= 

Defense  Mapping  Agency 

DMAAC 

Defense  Mapping  Agency  Aerospace  Center 

DMAHTC 

- 

Defense  Mapping  Agency  Hydrographic  Topographic  Center 

DoD 

= 

Department  of  Defense 

ECEF 

IB 

Earth-Centered,  Earth-Fixed 

ECI 

= 

Earth-Centered  Inertial 

ECM 

as 

Earth's  Center  of  Mass 

ED  50 

= 

European  Datum  1950 

ED  79 

European  Datum  1979 

EGM 

au 

Earth  Gravitational  Model 

FRG 

BK 

Federal  Republic  of  Germany 

GLONASS 

= 

Global  Navigation  Satellite  System 

GPS 

Global  Positioning  System 

GRS  80 

= 

Geodetic  Reference  System  1980 

lAG 

= 

International  Association  of  Geodesy 

lAU 

International  Astronomical  Union 

ITS 

ES 

Instantaneous  Terrestrial  System 

lUGG 

= 

International  Union  of  Geodesy  and  Geophysics 

MC&G 

SS 

Mapping,  Charting,  and  Geodesy 

MREs 

= 

Multiple  Regression  Equations 

NAD  27 

= 

North  American  Datura  1927 

NAD  83 

= 

North  American  Datura  1983 
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APPENDIX  E  (Cont'd) 
-ACRONYMS- 


NAVOCEANO 

Navstar  GPS 

NGS 

NNSS 

NSWC 

OSGB  36 
OSGB  SN  80 

PSAD  56 
RMS 

SAD  69 

TD 

TR 

USNO 

UK 

US 

UT 

VLB  I 
WGS 

WGS  72 
WGS  84 


Naval  Oceanographic  Office 
Navstar  Global  Positioning  System 
National  Geodetic  Survey 
Navy  Navigation  Satellite  System 

Naval  Surface  Warfare  Center  (formerly  Naval  Surface 
Weapons  Center) 

Ordnance  Survey  of  Great  Britain  1936 

Ordnance  Survey  of  Great  Britain  Scientific  Network 

1980 

Provisional  South  American  Datum  1956 

Root -Mean-Square 

South  American  Datiam  1969 

Tokyo  Datum 

Technical  Report 

United  States  Naval  Observatory 

United  Kingdom 

United  States 

Universal  Time 

Very  Long  Baseline  Interferometry 
World  Geodetic  System 
World  Geodetic  System  1972 
World  Geodetic  System  1984 
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